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Fifty years of developments in nuclear magnetic resonance (NMR) have resulted in an unrivaled 
degree of control of the dynamics of coupled two-level quantum systems. This coherent control of 
nuclear spin dynamics has recently been taken to a new level, motivated by the interest in quantum 
information processing. NMR has been the workhorse for the experimental implementation of 
quantum protocols, allowing exquisite control of systems up to seven qubits in size. Here, we 
survey and summarize a broad variety of pulse control and tomographic techniques which have 
been developed for and used in NMR quantum computation. Many of these will be useful in other 
quantum systems now being considered for implementation of quantum information processing 
tasks. 
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I. INTRODUCTION 

Precise and complete control of multiple coupled quan- 
tum systems is expected to lead to profound insights in 
physics as well as to novel applications, such as quantum 
computation (flsBlil)- Such coherent control is a major 
goal in atomic physics S IS 0)) quantum optics (0;13) 
and condensed matter research (Jy; |3 m LUj) ? but sur- 
prisingly, many of the leading experimental results are 
coming from one of the oldest areas of quantum physics: 
nuclear magnetic resonance (NMR). 

The development of NMR control techniques origi- 
nated in a strong demand for precise spectroscopy of 
complex molecules: NMR is the premier tool for pro- 
tein structure determination, and in modern NMR spec- 
troscopy, often thousands of precisely sequenced and 
phase controlled pulses are applied to molecules con- 
taining hundreds of nuclear spins. More recently, over 
the past seven years, a wide variety of complex quan- 
tum information processing tasks have been realized us- 
ing NMR, on systems ranging from two to seven quan- 
tum bits (qubits) in size, on molecules in liquid fill IT3 
Hll llll: llallrjl) . liquid crystal ifl^) . and solid state sam- 
ples j|l8t Il9j) . These demonstrations have been made 
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possible by application of a menagerie of new and pre- 
viously existing control techniques, such as simultaneous 
and shaped pulses, composite pulses, refocusing schemes, 
and effective Hamiltonians. These allow control and com- 
pensation for a variety of imperfections and experimen- 
tal artifacts invariably present in real physical systems, 
such as pulse imperfections, Bloch-Siegert shifts, unde- 
sired multiple-spin couplings, field inhomogeneities, and 
imprecise system Hamiltonians. 

The problem of control of multiple coupled quantum 
systems is a signature topic for NMR, and can be sum- 
marized as follows: given a system with Hamiltonian 
H = H sys + Hcontrolj where 7i sys is the Hamiltonian in 
the absence of any active control, and H con troi describes 
terms which are under external control, how can a de- 
sired unitary transformation U be implemented, in the 
presence of imperfections, and using minimal resources? 
Similar to other scenarios in which quantum control is a 
well-developed idea, such as in laser excitation of chemi- 
cal reactions ()20fl , 7i C ontroi arises from precisely timed se- 
quences of multiple pulses of electromagnetic radiation, 
applied phase-coherently, with different pulse widths, fre- 
quencies, phases, and amplitudes. However, importantly, 
in contrast to other areas of quantum control, in NMR 
H. sys is composed from multiple distinct physical pieces, 
i.e. the individual nuclear spins, providing the tensor 
product Hilbert space structure vital to quantum com- 
putation. Furthermore, the NMR systems employed in 
quantum computation are more well approximated as be- 
ing closed, as opposed to open, quantum systems. 

Nuclear spins and NMR provide a wonderful model and 
inspiration for the advance of coherent control over other 
coupled quantum systems, as many of the challenges and 
solutions are similar across the world of ato mic, molecu- 
lar, optical, and solid-state systems (see e.g. i|2lf) '). Here, 
we review the control techniques employed in the field 
of NMR quantum computation, focusing on methods 
which are robust under experimental implementation, 
and including experimental prescriptions for evaluation 
of the efficacy of the techniques. In contrast to other re- 
views |H H3) of and introductions HI El El) 
to NMR quantum computation which have appeared in 
the literature, we do not assume prior knowledge of, or 
give specialized descriptions of quantum computation al- 
gorithms, nor do we review NMR quantum computing 
experiments. And although we do not assume prior de- 
tailed knowledge of NMR, a self-contained treatment of 
several advanced topics, such as composite pulses, and 
refocusing, is included. Finally, as a primary purpose of 
this article is to elucidate control techniques which may 
generalize beyond NMR, we also assume a regime of op- 
eration in which relaxation and decoherence mechanisms 
are simple to treat and physical evolution is dominated 
by closed systems dynamics. 

The organization of this article is as follows. In sec- 
tion [HJ we briefly review the physics of NMR, using a 
Hamiltonian description of single and interacting nuclear 
spins-1/2 placed in a static magnetic field, controlled by 



radio-frequency fields. This establishes a foundation for 
the first major part of this review, section UTTl which dis- 
cusses the knobs the control Hamiltonian provides to con- 
struct all the elementary quantum gates, and the limita- 
tions that arise from the given system and control Hamil- 
tonian, as well as from instrumental imperfections. The 
second major part of this review, section IIVI presents 
three classes of advanced techniques for tailoring the con- 
trol Hamiltonian, which permit accurate quantum con- 
trol despite the existing limitations: the methods of am- 
plitude and frequency shaped pulses, composite pulses 
and average Hamiltonian theory. Finally, in section^ we 
conclude by describing a set of standard experiments, de- 
rived from quantum computation, which demonstrate co- 
herent qubit-control and can be used to characterize de- 
coherence. These include procedures for quantum state 
and process tomography, as well as methods to evaluate 
the fidelity of quantum states and gates. 

For further reading on NMR, we recommend the 
textbooks by Abragam (I29 r) . Ernst. Bodenhausen and 
Wokaun 1)301) and Slichter Ij31|) for rigorous discussions 
of the nuclear spin Hamiltonian and standard pulse se- 
quences, Freeman ll.32l) for an intuitive explanation of ad- 
vanced techniques for control of the spin evolution, and 
Levitt (33) for an intuitive understanding of the physics 
underlying the spin dynamics. Many useful reviews on 
specific NMR techniques are compiled in the Encyclope- 
dia of NMR Q . 

For additional reading on quantum computation, we 
recommend the book by Nielsen and Chuang |2j) for the 
basic theory of quantum information and computation 
and Refs. Q l35h for a review of the state of the art in 
experimental quantum information processing. Ref. l)36|) 
gives a simple introduction to quantum computation. 
Excellent presentations of quantum algorithms are given 
in Refs. JHIEI)- 

The original papers introducing NMR quantum com- 
puting are Refs (IM III El EI Refs. JH EI give 
elementary introductions to NMR quantum computing, 
and introductions geared towards NMR spectroscopists 
are presented in Refs l|25t l2t?) . Summaries of NMR quan- 
tum c omp uting experiments and techniques are given in 

Refs. mM m>- 



II. THE NMR SYSTEM 

We begin with a description of the NMR system, based 
on its system Hamiltonian and the control Hamiltonian. 
The system Hamiltonian gives the energy of single and 
coupled spins in a static magnetic field, and the con- 
trol Hamiltonian arises from the application of radio- 
frequency pulses to the system at, or near, its resonant 
frequencies. A rotating reference frame is employed, pro- 
viding a very convenient description. 
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A. The system Hamiltonian 

1. Single spins 




The time evolution of a spin- 1/2 particle (we will not 
consider higher order spins in this paper) in a magnetic 
field Bq along z is governed by the Hamiltonian 



Ho = -h-yB I z = -hu I z = 



~hhj /2 
huj a /2 



, (1) 

where 7 is the gyromagnetic ratio of the nucleus, loq/2tt is 
the Larmor frequency 1 and I z is the angular momentum 
operator in the z direction. I z , I x , and I y relate to the 
well-known Pauli matrices as 



a x = 21, 
where, in matrix notation, 



(Jy — 2I y , 
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The interpretation of Eq.^is that the |0) or |f) energy 
(given by (0|W|0), the upper left element of 7i) is lower 
than the |1) or ||) energy by an amount huj n , as 

illustrated in the energy diagram of Fig. ^ The energy 
splitting is known as the Zccman splitting. 
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FIG. 2 Precession of a spin- 1/2 about the axis of a static 
magnetic field. 



The nuclear spin Hamiltonian for a molecule with n 
uncoupled nuclei with is thus given by 



Wo = - - fiihiBo ii = -J2 huj o n 



where the i superscripts label the nuclei. 



(4) 



nucleus 


X H 2 H 13 C 


15 N 


19p 31p 


UlQ /2-7T 


500 77 126 


-51 


470 202 



TABLE I Larmor frequencies [MHz] for some relevant nuclei, 
at 11.74 Tesla. 




FIG. 1 Energy diagram for a single spin-1/2. 

We can pictorially understand the time evolution U = 
e -iHt/h un( j er th. e Hamiltonian of Eq. ^ as a process- 
ing motion of the Bloch vector about Bq, as shown in 
Fig. [2] As is conventional, we define the z axis of the 
Bloch sphere as the quantization axis of the Hamiltonian, 
with |0) along +z and |1) along —z. 

For the case of liquid-state NMR, which we will largely 
restrict ourselves to in this article, typical values of B 
are 5—15 Tesla, resulting in precession frequencies u>q of 
a few hundred MHz, the radio-frequency range. 

Spins of different nuclear species (heteronuclear spins) 
can be easily distinguished spectrally, as they have very 
distinct values of 7 and thus also very different Larmor 
frequencies (Table QJ. Spins of the same nuclear species 
(homonuclear spins) which are part of the same molecule 
can also have distinct frequencies, by amounts known as 
their chemical shifts, di. 



1 We will sometimes leave the factor of 2-k implicit and call ljq the 
Larmor frequency. 




FIG. 3 Fluorine NMR spectrum (absolute value) centered 
around ~ 470 MHz of a specially designed molecule (shown 
on the right). The five main lines in the spectrum correspond 
to the five fluorine nuclei in the molecule. The two small lines 
derive from impurities in the sample. The NMR spectra were 
acquired by recording the oscillating magnetic field produced 
by a large ensemble of precessing spins, and by taking the 
Fourier transform of this time-domain signal. The precession 
motion of the spins is started by applying a radio-frequency 
pulse (section IlLB.H which tips the spins from their equilib- 
rium position along the z axis into the x — y plane. 

The chemical shifts arise from partial shielding of the 
externally applied magnetic field by the electron cloud 
surrounding the nuclei. The amount of shielding depends 
on the electronic environment of each nucleus, so like nu- 
clei with inequivalent electronic environments have dif- 
ferent chemical shifts. Pronounced asymmetries in the 
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molecular structure generally promote strong chemical 
shifts. The range of typical chemical shifts (7j varies from 
nucleus to nucleus, e.g. f=a 10 parts per million (ppm) 
for 1 H, « 200 ppm for 19 F and « 200 ppm for 13 C. At 
Bq = 10 Tesla, this corresponds to a few kHz to tens of 
kHz (compared to cJn's of several hundred MHz). As an 
example, Fig. [3] shows an experimentally measured spec- 
trum of a molecule containing five fluorine spins with 
inequivalent chemical environments. 

In general, the chemical shift can be spatially 
anisotropic and must be described by a tensor. In liquid 
solution, this anisotropy averages out due to rapid tum- 
bling of the molecules. In solids, the anisotropy means 
that the chemical shifts depend on the orientation of the 
molecule with respect to Bq. 

2. Interacting spins 

For nuclear spins in molecules, nature provides two dis- 
tinct interaction mechanisms which we now describe, the 
direct dipole-dipole interaction, and the electron medi- 
ated Fermi contract interaction known as J-coupling. 

Direct coupling. The magnetic dipole-dipole interac- 
tion is similar to the interaction between two bar mag- 
nets in each other's vicinity. It takes place purely through 
space — no medium is required for this interaction — and 
depends on the internuclear vector fy connecting the two 
nuclei i and j, as described by the Hamiltonian 

(5) 

where /in is the usual magnetic permeability of free space 
and P is the magnetic moment vector of spin i. This 
expression can be progressively simplified as various con- 
ditions are met. These simplifications rest on averaging 
effects and can be explained within the general frame- 
work of average-Hamiltonian theory fsection llV.C(l . 

For large ujq = j l B (i.e. at high B ), Hd can be 
approximated as 

n ° = E fSSa - 3cos2 %)t 3/ ^' - r ■ fj ] . (e) 

i<j ' ^ ' 

where % is the angle between B and r^. When |w — uJq\ 
is much larger than the coupling strength, the transverse 
coupling terms can be dropped, soHd simplifies further 
to 

W^ESVs-rf^, (7) 

i<j ' ' 

which has the same form as the J-coupling we describe 
next (Eq.EJl. 

For molecules in liquid solution, both intramolecular 
dipolar couplings (between spins in the same molecule) 
and intermolecular dipolar couplings (between spins in 



different molecules) are averaged away due to rapid tum- 
bling. This is the case we shall focus on in this article. In 
solids, similarly simple Hamiltonians can be obtained by 
applying multiple-pulse sequences which average out un- 
desired coupling terms l)43T) . or by physically spinning the 
sample at an angle of arccos(l/v / 3) (the "magic angle") 
with respect to the magnetic field. 

Indirect coupling. The second interaction mech- 
anism between nuclear spins in a molecule is the J- 
coupling or scalar coupling. This interaction is mediated 
by the electrons shared in the chemical bonds between 
the atoms, and due to the overlap of the shared electron 
wavefunction with the two coupled nuclei, a Fermi con- 
tact interaction. The through-bond coupling strength J 
depends on the respective nuclear species and decreases 
with the number of chemical bonds separating the nuclei. 
Typical values for J are up to a few hundred Hz for one- 
bond couplings and down to only a few Hertz for three- 
or four-bond couplings. The Hamiltonian is 

Hj = hJ22*Jiji l -i j - /^2t./, ; :/;/; • /;/• ■ 

i<j i<j 

where is the coupling strength between spins i and j. 
Similar to the case of dipolar coupling, Eq. [8] simplifies 
to 

n 

n. ] = hY J 2^Wi, (9) 

when \u>i — uij\ 2ir\Jij\, a condition easily satisfied for 
heteronuclear spins and which can also be satisfied for 
small homonuclear molecules. 

The interpretation of the scalar coupling term of Eq. [5] 
is that a spin "feels" a static magnetic field along ±z 
produced by the neighboring spins, in addition to the 
externally applied Bq field. This additional field shifts 
the energy levels as in Fig. As a result, the Larmor 
frequency of spin i shifts by — Jy/2 if spin j is in |0) and 
by +Jij/2 if spin j is in |1). 




FIG. 4 Energy level diagram for (dashed lines) two uncoupled 
spins and (solid lines) two spins coupled by a Hamiltonian of 
the form of Eq. Qor Eq.|5](in units of K). 

In a system of two coupled spins, the frequency spec- 
trum of spin i therefore actually consists of two lines 
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separated by Jy and centered around u)q, each of which for the static field Bo, 
can be associated with the state of spin j, |0) or For 
three pairwise coupled spins, the spectrum of each spin 
contains four lines. For every additional spin, the number 
of lines per multiplet doubles, provided all the couplings 
are resolved and different lines do not lie on top of each 
other. This is illustrated for a five spin system in Fig. [S] 
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FIG. 5 The spectrum of spin Fi in the molecule of Fig. El 
This is an expanded view of the left line in the spectrum of 
Fig. El Frequencies are given with respect to Uq. The state 
of the remaining spins is as indicated, based on J12 < and 
J13, J14, J15 > 0; furthermore, |Ji 2 | > IJ13I > |Jis| > |Ji4|- 



The magnitude of all the pairwise couplings can be 
found by looking for common splittings in the multi- 
plets of different spins. The relative signs of the J cou- 
plings can be determined via appropriate spin-selective 
two-pulse sequences, known in NMR as two-dimensional 
correlation (soft-COSY) experiments ifrU) or via line- 
selective continuous irradition; both approaches are re- 
lated to the cnot gate ( section 1111. A. 3 [I . The signs can- 
not be obtained from just the simple spectra. 

In summary, the simplest form of the Hamiltonian for 
a system of n coupled nuclear spins is thus (from Eqs. 0] 
and© 



H sys = - £ K<4 II + h^^Wi 



(10) 



l<3 



In almost all NMR quantum computing experiments per- 
formed to date, the system is well described by a Hamil- 
tonian of this form. 



B. The control Hamiltonian 

1. Radio-frequency fields 

We turn now to physical mechanisms for controlling 
the NMR system. The state of a spin- 1/2 particle in a 
static magnetic field Bq along z can be manipulated by 
applying an electromagnetic field B\(t) which rotates in 
the x — y plane at uj r f, at or near the spin precession fre- 
quency loq. The single-spin Hamiltonian corresponding 
to the radio-frequency (RF) field is, analogous to Eq. 2] 



H r f = —tvyBi cos(u> r ft + 4>)I X + sin(u> r ft + 4>)I y 

'(11) 

where </> is the phase of the RF field, and B\ its ampli- 
tude. Typical values for u>\ = 7B1 are up to ~ 50 kHz in 
liquid NMR and up to a few hundred kHz in solid NMR 
experiments. For n spins, we have 



H r f = — V] fcyiBi cos(uj r ft + <j>)I x + sm(oj r ft + </>)/* 

In practice, a magnetic field is applied which oscillates 
along a fixed axis in the laboratory, perpendicular to the 
static magnetic field. This oscillating field can be de- 
composed into two counter-rotating fields, one of which 
rotates at u> r f in the same direction as the spin and so 
can be set on or near resonance with the spin. The other 
component rotates in the opposite direction and is thus 
very far off-resonance (by about 2cJn). As we shall see, its 
only effect is a negligible shift in the Larmor frequency, 
called the Bloch-Siegert shift lf45T) . 

Note that both the amplitude B\ and phase <f> of the 
RF field can be varied with time 2 , unlike the Larmor 
precession and the coupling terms. As we will shortly 
see, it is the control of the RF field phases, amplitudes, 
and frequencies, which lie at the heart of quantum control 
of NMR systems. 



2. The rotating frame 

The motion of a single nuclear spin subject to both 
a static and a rotating magnetic field is rather complex 
when described in the usual laboratory coordinate sys- 
tem (the lab frame). It is much simplified, however, by 
describing the motion in a coordinate system rotating 
about z at to r f (the rotating frame) : 



\^) rot = exp(-iuj r ftl z )\ip) . 



(13) 



dm 



Substitution of \ip) in the Schrodinger equation ifr^jj 
H\ib) with 



TL = —hujQ l z — fujj\ 



cos(w r jt + <j))I x + $m(uj r ft - 



gives — = TL rot \ip) ro , where 

H rot — —h(ujQ—uj r f) I z — huJi cos<f)I x -\-sm.(f)Iy 



(14) 



(15) 



2 For example, the Varian Instruments Unity Inova 500 NMR spec- 
trometer achieves a phase resolution of 0.5° and has 4095 linear 
steps of amplitude control, with a time-base of 50 ns. Additional 
attenuation of the amplitude can be done on a logarithmic scale 
over a range of about 80 dB, albeit with a slower timcbasc. 
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Naturally, the RF field lies along a fixed axis in the frame 
rotating at tu r f. Furthermore, if Lu r f — lu , the first term 
in Eq. If 51 vanishes. In this case, an observer in the ro- 
tating frame will see the spin simply precess about B\ 
(Fig. Ek) , a motion called nutation. The choice of <j> con- 
trols the nutation axis. An observer in the lab frame sees 
the spin spiral down over the surface of the Bloch sphere 
(Fig.Ep). 





(a) 



FIG. 6 Nutation of a spin subject to a transverse RF field (a) 
observed in the rotating frame and (b) observed in the lab 
frame. 

If the RF field is off-resonance with respect to the spin 
frequency by Ao; = luq — u) r f, the spin precesses in the 
rotating frame about an axis tilted away from the z axis 
by an angle 



a — arctan(u;i/Au/) , 
and with frequency 



as illustrated in Fig. 



(16) 



(17) 




FIG. 7 Axis of rotation (in the rotating frame) during an 
off-resonant radio-frequency pulse. 

It follows that the RF field has virtually no effect on 
spins which are far off-resonance, since a is very small 
when |Aw| ^> u>i (see Fig. [SJ. If all spins have well- 
separated Larmor frequencies, we can thus in principle se- 
lectively rotate any one qubit without rotating the other 
spins. 

Moderately off-resonance pulses (|Aw| w u>i) do rotate 
the spin, but due to the tilted rotation axis, a single such 
pulse cannot, for instance, flip a spin from |0) to |1) (see 
again Fig.|HJ). Of course, off-resonance pulses can also be 
useful, for instance for direct implementation of rotations 
about an axis outside the x — y plane. 



FIG. 8 Trajectory in the Bloch sphere described by a qubit 
initially in |0) (along +z), after a 250 us pulse of strength 
Mi = 1 kHz is applied off-resonance by 0, 0.5, 1, ...4 kHz. 
On-resonance, the pulse produces a 90° rotation. Far off- 
resonance, the qubit is hardly rotated away from |0). 



We could also choose to work in a frame rotating at u>q 
(instead of 0J r f), where 



-fkoi 



cos((w r/ - U) )t + (f>)I x 



+ SUx((oJ r f — OJo)t + 4>)I y 



(18) 



This transformation does not give a convenient time- 
independent RF Hamiltonian (unless Lu r f — ujq), as was 
the case for 7i rot in Eq.^J However, it is a natural start- 
ing point for the extension to the case of multiple spins, 
where a separate rotating frame can be introduced for 
each spin: 



Uexp(-iu;5t/*) 



(19) 



In the presence of multiple RF fields indexed r, the RF 
Hamiltonian in this multiply rotating frame is 



n rot - ^-M[cos(K / -c4)i + ^)4 



(20) 



where the amplitudes uj\ and phases <jf are under user 
control. 

The system Hamiltonian of Eg. HOI is simplified, in the 
rotating frame of Eq. ^3 the 1\ terms drop out leaving 
just the JijP z P z couplings, which remain invariant. Note 
that coupling terms of the form P ■ I 3 do not transform 
cleanly under Eq. 1191 

Summarizing, in the multiply rotating frame, the NMR 



Hamiltonian Tt = 7Y S 



control 



takes the form 



(21) 



^control = ^-hhj{ CO$,((uj r rf ~ Ujl)t + (j) r )P x 



+ Bin((w£-4)t + ^)j; 



(22) 



7 



C. Relaxation and decoherence 



III. ELEMENTARY PULSE TECHNIQUES 



One of the strengths of nuclear spins as quantum bits 
is precisely the fact that the system is very well isolated 
from the environment, allowing coherence times to be 
long compared with the dynamical timescales of the sys- 
tem. Thus, our discussion here focuses on closed system 
dynamics, and it is important to be aware of the limits 
of this approximation. 

The coupling of the NMR system to the environment 
may be described by an additional Hamiltonian term 
Ticnv, whose magnitude is small compared to that of 7i sys 
or Hcontroi- It is this coupling which leads to decoherence, 
the loss of quantum information, which is traditionally 
parameterized by two rates: T\, the energy relaxation 
rate, and T2, the phase randomization rate (see also Sec- 
tions IV.A.4I and IV.A.5)) . 

T2 originates from spin-spin couplings which are imper- 
fectly averaged away, or unaccounted for in the system 
Hamiltonian. For example, in molecules in liquid solu- 
tion, spins on one molecule may have a long range, weak 
interaction with spins on another molecule. Fluctuat- 
ing magnetic fields, caused by spatial anisotropy of the 
chemical shift, local paramagnetic ions, or unstable labo- 
ratory fields, also contribute to T^. Nevertheless, in well 
prepared samples and in a good experimental apparatus 
at reasonably high magnetic fields, the T2 for molecules 
in solution is easily on the order of one second or more. 
This decoherence mechanism can be identified with elas- 
tic scattering in other physical systems; it does not lead 
to loss of energy from the system. 

T\ originates from couplings between the spins and the 
"lattice," that is, excitation modes which can carry away 
energy quanta on the scale of the Larmor frequency. For 
example, these may be vibrational quanta, paramagnetic 
ions, chemical reactions such as ions exchanging with the 
solvent, or spins with higher order magnetic moments 
(such as 2 H, 17 C1, or 35 Br) which relax quickly due to 
their quadrupolar moments interacting with electric field 
gradients. In well chosen molecules and liquid samples 
with good solvents, T\ can easily be tens of seconds, while 
isolated nuclei embedded in solid samples with a spin- 
zero host crystal matrix (such as 31 P in 28 Si) can have T\ 
times of days. This mechanism is analogous to inelastic 
scattering in other physical systems. 

The description of relaxation in terms of only two pa- 
rameters is known to be an oversimplification of reality, 
particularly for coupled spin systems, in which coupled 
relaxation mechanisms appear (|4a : |47j). Nevertheless, 
the independent spin decoherence model is useful for its 
simplicity and because it can capture well the main ef- 
fects of decoherence on simple NMR quantum compu- 
tations fllrih . which are typically designed as pulse se- 
quences shorter in time than T^. 



This section begins our discussion of the main subject 
of this article, a review of the control techniques devel- 
oped in NMR quantum computation for coupled two- 
level quantum systems. We begin with a quick overview 
of the language of quantum circuits and its important 
universality theorems, then connect this with the lan- 
guage of pulse sequences as used in NMR, and indicate 
how pulse sequences can be simplified. The main approx- 
imations employed in this section are that pulses can be 
strong compared with the system Hamiltonian while se- 
lectively addressing only one qubit at a time, and can be 
perfectly implemented. The limits of these approxima- 
tions are discussed in the last part of the section. 



A. Quantum control, quantum circuits, and pulses 

The goal of quantum control, in the context of quan- 
tum computation, is the implementation of a unitary 
transformation U, specified in terms of a sequence U — 
UkUk-i ■ ■ ■ U2U1 of standard "quantum gates" Ui, which 
act locally (usually on one or two qubits) and are simple 
to implement. As is conventional for unitary operations, 
the Ui are ordered in time from right to left. 



1. Quantum gates and circuits 

The basic single-qubit quantum gates are rotations, 
defined as 



Rn{0) = exp 



i9n ■ a 



(23) 



where n is a (three-dimensional) vector specifying the 
axis of the rotation, is the angle of rotation, and 
3 = a x x + a y y + tr z z is a vector of Pauli matrices. It 
is also convenient to define the Pauli matrices (see Eq. |3J 
themselves as logic gates, in terms of which a x can be 
understood as being analogous to the classical not gate, 
which flips |0) to |1) and vice versa. In addition, the 
Hadamard gate H and 7r/8 gate T 



H 



1 

71 



1 1 


, T = 


1 -1 





1 

exp(i7r/4) 



(24) 



are useful and widely employed. These, and any other 
single qubit transformation U can be realized using a 
sequence of rotations about just two axes, according to 
Bloch's theorem: for any single-qubit U, there exist real 
numbers a, (3, 7 and 5 such that 

U = e la R x ((3)R y (j)R x (S). (25) 

The basic two-qubit quantum gate is a controlled-NOT 
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(cnot) gate 



C^cnot — 



10 
10 
1 
10 



(26) 



where the basis elements in this notation are 1 00) , |01), 
1 10), and 1 11) from left to right and top to bottom. 
t^CNOT flips the second qubit (the target) if and only 
if the first qubit (the control) is |1). This gate is 
the analogue of the classical exclusive-OR gate, since 
Ucnot\x,u} = \x,x®y), for x,y € {0,1} and where 
© denotes addition modulo two. 

A basic theorem of quantum computation is that up to 
an irrelevant overall phase, any U acting on n qubits can 
be composed from ZTcnot and Rh(6) gates (@). Thus, 
the problem of quantum control can be reduced to im- 
plementing [/cnot an d single qubit rotations, where at 
least two non-trivial rotations are required. Other such 
sets of universal gates are known, but this is the one 
which has been employed in NMR. 

These gates and sequences of such gates may be con- 
veniently represented using quantum circuit diagrams, 
employing standard symbols. We shall use a notation 
commonly employed in the literature (|2j) in this article. 



2. Implementation of single qubit gates 

Rotations on single qubits may be implemented di- 
rectly in the rotating frame using RF pulses. From the 
control Hamiltonian, Eq. [23 it follows that when an RF 
field of amplitude u>i is applied to a single-spin system at 
LU r f = w 0j t ne spin evolves under the transformation 



U = exp iuji (cos 4>I X + sin (f> I y )t pw 



(27) 



where t pw is the pulse width (or pulse length), the time 
duration of the RF pulse. U describes a rotation in the 
Bloch sphere over an angle 9 proportional to the product 
of t pw and uj\ = jB±, and about an axis in the x — y 
plane determined by the phase <f>. 

Thus, a pulse with phase <f> — 71 an d oj\t pw — it/2 
will perform R x (90) (see Ea. l23[) . which is a 90° rotation 
about x, denoted for short as X. A similar pulse but 
twice as long realizes a R x (180) rotation, written for short 
as X 2 . By changing the phase of the RF pulse to tf> — 
— 7r/2, Y and Y 2 pulses can similarly be implemented. 
For tj> — 0, a negative rotation about x, denoted R x {— 90) 
or X, is obtained, and similarly <j> — it/2 gives Y. For 
multi-qubit systems, subscripts are used to indicate on 
which qubit the operation acts, e.g. Z\ is a 180° rotation 
of qubit 3 about — z. 

It is thus not necessary to apply the RF field along 
different spatial axis in the lab frame to perform x and y 
rotations. Rather, the phase of the RF field determines 
the nutation axis in the rotating frame. Furthermore, 



note that only the relative phase between pulses applied 
to the same spin matters. The absolute phase of the first 
pulse on any given spin does not matter in itself. It just 
establishes a phase reference against which the phases of 
all subsequent pulses on that same spin, as well as the 
read-out of that spin, should be compared. 

We noted earlier that the ability to implement arbi- 
trary rotations about x and y is sufficient for performing 
arbitrary single-qubit rotations (Eq. I25fl . Since z rota- 
tions are very common, two useful explicit decomposi- 
tions of R z {9) in terms of x and y rotations are: 



R z {6) = XRy(0)X = YR x (-0)Y . 



3. Implementation of two-qubit gates 



(28) 



The most natural two-qubit gate is the one generated 
directly by the spin-spin coupling Hamiltonian. For nu- 
clear spins in a molecule in liquid solution, the coupling 
Hamiltonian is given by Eq. [5] (in the lab frame as well 
as in the rotating frame) , from which we obtain the time 
evolution operator Uj(t) = exp[—i2irJIgI^t], or in ma- 
trix form 



Uj(t) 



-inJt/2 o 

e +if.«/2 o 
e +if.«/2 o 

e -i*Jt/z 

(29) 

Allowing this evolution to occur for time t = 1/2 J gives 
a transformation known as the controlled phase gate, up 
to a 90° phase shift on each qubit and an overall (and 
thus irrelevant) phase: 



C^CPHASE = V^lZ\Z 2 Uj{l/2J) 



10 
10 
10 
-1 



(30) 

This gate is equivalent to the well-known CNOT gate up 
to a basis change of the target qubit and a phase shift on 
the control qubit: 



C^CNOT 



iZjY 2 



iZ 1 Z 2 Uj{l/2J) 



Yo 



ViZ 1 Z 2 X 2 Uj(l/2J)Y 2 

1 0" 

10 

1 

10 



(31) 



The core of this sequence, X 2 Uj(l/2 J)Y 2 , can be graph- 
ically understood via Fig. ED (4fJ), assuming the spins 
start along ±z. First, a spin-selective pulse on spin 2 
about y (an rf pulse centered at Wq/27t and of a spec- 
tral bandwidth such that it covers the frequency range 
u!q/2tt±Ji 2 /2 but not 0Jq/2it± Ji 2 /2, rotates spin 2 from 
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FIG. 9 Bloch-sphere representation of the operation of the 
CNOT12 gate between two qubits 1 and 2 coupled by 
h2irJlllz. Here, qubit 2 starts off in |0) (along z) and is de- 
picted in a reference frame rotating about z at loq/2it. Solid 
and dashed arrows correspond to the case where qubit 1 is |0) 
and |1) respectively. 



z to x. Next, the spin system is allowed to freely evolve 
for a duration of 1/2 Ji 2 seconds. Because the precession 
frequency of spin 2 is shifted by ±J 12 /2 depending on 
whether spin 1 is in |1) or |0) (see Fig. QJ, spin 2 will 
arrive in 1/2 J seconds at either +y or — y, depending on 
the state of spin 1. Finally, a 90° pulse on spin 2 about 
the x axis rotates spin 2 back to +z if spin 1 is |0), or to 
— z if spin 1 is in |1). 

The net result is that spin 2 is flipped if and only if 
spin 1 is in |1), which corresponds exactly to the classical 
truth table for the CNOT. The extra z rotations in Eg. 1311 
are needed to give all elements in (7cnot the same phase, 
so the sequence works also for superposition input states. 

An alternative implementation of the CNOT gate, up 
to a relative phase factor, consists of applying a line- 
selective 180° pulse at ujq + J12/2 (see Fig.0J. This pulse 
inverts spin 2 (th e target qubit) if and only if spin 1 (the 
control) is 1 1) l|42|) . In general, if a spin is coupled to more 
than one other spin, half the lines in the multiplet must 
be selectively inverted in order to realize a CNOT. Ex- 
tensions to doubly-controlled nots are straightforward: 
in a three-qubit system for example, this can be realized 
through inversion of one out of the eight lines |48lh As 
long as all the lines are resolved, it is in principle possible 
to invert any subset of the lines. Demonstrations using 
very long multi-frequency pulses have been performed 
with up to five qubits (|49j). However, this approach can- 
not be used whenever the relevant lines in the multiplet 
fall on top of each other. 

If the spin-spin interaction Hamiltonian is not of the 
form P Z P Z but contains also transverse components (as in 
Eqs.[SHS]and[SJ), other sequences of pulses are needed to 
perform the CPHASE and CNOT gates. These sequences 
are somewhat more complicated l)50fl . 

If two spins are not directly coupled two each other, it 
is still possible to perform a CNOT gate between them, as 
long as there exists a network of couplings that connects 
the two qubits. For example, suppose we want to perform 
a CNOT gate with qubit 1 as the control and qubit 3 
as the target, CNOT13, but 1 and 3 are not coupled to 
each other. If both are coupled to qubit 2, as in the 
coupling network of Fig. ^] (b), we can first swap the 
state of qubits 1 and 2 (via the sequence CNOT12 CNOT21 
CNOT12), then perform a CNOT23, and finally swap qubits 
1 and 2 again (or relabel the qubits without swapping 



back). The net effect is CNOT13. By extension, at most 
0(n) swap operations are required to perform a CNOT 
between any pair of qubits in a chain of n spins with just 
nearest- neighbor couplings (Fig. Hub ) . swap operations 
can also be used to perform two-qubit gates between any 
two qubits which are coupled to a common "bus" qubit 

(Fig. not). 




(c) 



FIG. 10 Three possible coupling networks between five 
qubits. (a) A full coupling network. Such networks will in 
practice always be limited in size, as physical interactions tend 
to decrease with distance, (b) A nearest-neighbor coupling 
network. Such linear chains with nearest-neighbor couplings, 
or two-dimensional variants are used in many solid-state pro- 
posals, (c) Coupling via a "bus" . This is the case of ion-trap 
schemes for example. Similar to case (a), the bus degree of 
freedom will in reality couple well to only a finite number of 
qubits. 

Conversely, if a qubit is coupled to many other qubits 
(Fig. Hub ) and we want to perform a CNOT between just 
two of them, we must remove the effect of the remaining 
couplings. This can accomplished using the technique of 
refocusing, which has been widely adopted in a variety 
of NMR experiments. 

4. Refocusing: turning off undesired couplings 

The effect of coupling terms during a time interval 
of free evolution can be removed via so-called "refocus- 
ing" techniques. For coupling Hamiltonians of the form 
as is often the case in liquid NMR experiments (see 
Eq. , the refocusing mechanism can be understood at 
a very intuitive level. Reversal of the effect of coupling 
Hamiltonians of other forms, such as in Eqs. Ellwand 00 
is less intuitive, but can be understood within the frame- 
work of average Hamiltonian theory (section llV.C|) . 

Let us first look at two ways of undoing in a two- 
qubit system. In Fig. Illb . the evolution of qubit 1 in 
the first time interval r is reversed in the second time 
interval, due to the 180° pulse on qubit 2. In Fig. ITTb . 
qubit 1 continues to evolve in the same direction all the 
time, but the first 180° pulse causes the two components 
of qubit 1 to be refocused by the end of the second time 
interval. The second 180° pulse ensures that both qubits 
always return to their initial state. 

Mathematically, we can see how refocusing of J cou- 
plings works using the fact that for all r 

X\ Uj{t) Xl = Uj(-t) = X\ Uj(t) X\ , (32) 
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(a)^J°Ldelay(T 




FIG. 11 Bloch-sphere representation of the operation of two 
simple schemes to refocus the coupling between two coupled 
qubits. The diagram shows the evolution of qubit 1 (in the ro- 
tating frame) initially along —y, when qubit 2 is in |0) (solid) 
or in |1) (dashed). The refocusing pulse can be applied to 
either (a) qubit 2 or (b) qubit 1. 



which leads to 

X\ Uj (r) X\ U.j (r) = I = X\ Uj (r) X\ Uj (r) . (33) 

Replacing all Xf with Y?, the sequence works just the 
same. However, if we use sometimes X? and sometimes 
Y^ 2 , we get the identity matrix only up to some phase 
shifts. Also, if we applied pulses on both qubits simul- 
taneously, e.g. XiX%Uj(t) X\X^~U j(t), the coupling 
would not be removed. 

Fig.^Jgives insight in refocusing techniques in a multi- 
qubit system. Specifically, this scheme preserves the ef- 
fect of J12, while effectively inactivating all the other cou- 
plings. The underlying idea is that a coupling between 
spins i and j acts "forward" during intervals where both 
spins have the same sign in the diagram, and acts "in 
reverse" whenever the spins have opposite signs. When- 
ever a coupling acts forward and in reverse for the same 
duration, it has no net effect. 



1 



+ + + + 



+ + + 



2 
3 
4 

FIG. 12 Refocusing scheme for a four spin system, designed 
to leave J12 active the whole time but to neutralize the effect 
of the other Jij. The interval is divided into slices of equal 
duration, and the "+" and "— " signs indicate whether a spin 
is still in its original position, or upside down. The black 
rectangles represent 180° pulses, which flip the corresponding 
spin. 

Systematic methods for designing refocusing schemes 
for multi-qubit systems have been developed specifically 
for the purpose of quantum computing. The most com- 
pact scheme is based on Hadamard matrices l|5lt 15^) . 
A Hadamard matrix of order n, denoted by H(n), is an 
n x n matrix with entries ±1, such that 



The rows are thus pairwise orthogonal, and any two rows 
agree in exactly half of the entries. Identifying +1 and 
— 1 with + and — as in the diagram of Fig. ^1 we see 
that H(n) gives a valid decoupling scheme for n spins 
using only n time intervals. An example of H{\2) is 



+ 
+ 



+ 



+ 
+ 



1_ 

h h H 

H 1 1 

h h h + 

1 h + + 



h 

H 

1 

+ h 

1 

+ + + 

+ H 

h + 

h 



(35) 



If we want the coupling between one pair of qubits to 
remain active while removing the effect of all other cou- 
plings, we can simply use the same row of H(n) for those 
two qubits. 

H(n) does not exist for all n, but we can always find 
a decoupling sequence for n qubits by taking the first 
n rows of -ff(n), with n the smallest integer that sat- 
isfies n > n with known H(n). From the properties of 
Hadamard matrices, we can show that n/n is always close 
to 1 l)5 ijh So decoupling schemes for n spins require n 
time intervals and no more than nn 180° pulses. 

Another systematic approach to refocusing sequences 
is illustrated via the following 4-qubit scheme |53}) : 



+ + 
+ + 
+ + 



+ 



+ 



(36) 



H(n)H(n) 



nl . 



(34) 



For every additional qubit, the number of time intervals 
is doubled, and 180° pulses are applied to this qubit after 
the first, third, fifth, . . . time interval. The advantage of 
this scheme over schemes based on Hadamard matrices, 
is that it does not require simultaneous rotations of mul- 
tiple qubits. The main drawback is that the number of 
time intervals increases exponentially. 

We end this subsection with three additional remarks. 
First, each qubit will generally be coupled to no more 
than a fixed number of other qubits, since coupling 
strengths tend to decrease with distance. In this case, all 
refocusing schemes can be greatly simplified i|51r 153: l53f) . 

Second, if the forward and reverse evolutions under 
Jij are not equal in duration, a net coupled evolution 
takes place corresponding to the excess forward or reverse 
evolution. In principle, therefore, we can organize any 
refocusing scheme such that it incorporates any desired 
amount of coupled evolution for each pair of qubits. 



11 



Third, refocusing sequences can also be used to remove 
the effect of P z terms in the Hamiltonian. Of course, 
these terms vanish in principle if we work in the mul- 
tiply rotating frame (see Eg. 1211 . However, there may 
be some spread in the Larmor frequencies, for instance 
due to magnetic field inhomogeneities. This effect can 
then be reversed using refocusing pulses, as is routinely 
accomplished in spin-echo experiments (section IV.A.4fl . 



5. Pulse sequence simplification 

There are many possible pulse sequences which in an 
ideal world result in exactly the same unitary transfor- 
mation. Good pulse sequence design therefore attempts 
to find the shortest and most effective pulse sequence that 
implements the desired transformations. In section ITvl 
we will see that the use of more complex pulses or pulse 
sequences may sometimes increase the degree of quantum 
control. Here, we look at three levels of pulse sequence 
simplification. 

At the most abstract level of pulse sequence simplifi- 
cation, careful study of a quantum algorithm can give 
insight in how to reduce the resources needed. For ex- 
ample, a key step in both the modified Deutsch-Jozsa 
algorithm JHj) and the Grover algorithm |5f|) can be de- 
scribed as the transformation — > \x)\x © y), where 
\y) is set to (|0) — |l))/v2, so that the transformation in 
effect is |z)(|0) - |1»/V2 - (-r/W|z>(|0) - |1»/V2. 
We might thus as well leave out the last qubit as it is 
never changed. 

At the next level, that of quantum circuits, we can use 
simplification rules such as those illustrated in Fig. ED In 
this process, we can fully take advantage of commutation 
rules to move building blocks around, as illustrated in 
Fig. 1141 Furthermore, gates which commute with each 
other can be executed simultaneously. Finally, we can 
take advantage of the fact that most building blocks have 
many equivalent implementations, as shown for instance 
in Fig. E3 
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FIG. 14 Commutation of unitary operators can help simplify 
quantum circuits by moving building blocks around such that 
cancellation of operations as in Fig. 1131 become possible. For 
example, the three segments (separated by dashed lines) in 
these two equivalent realizations of the toffoli gate (doubly- 
controlled not) commute with each other and can thus be 
executed in any order. 
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FIG. 15 Choosing one of several equivalent implementations 
can help simplify quantum circuits, again by enabling can- 
cellation of operations as in Fig. 1131 For instance, the two 
control qubits in the toffoli gate play equivalent roles, so 
they can be interchanged. 



qubits may still be along ±i, in which case their mutual 
P Z I 3 Z couplings have no effect and thus need not be refo- 
cused. Similarly, if a subset of the qubits can be traced 
out at some point in the sequence, the mutual interac- 
tion between these qubits does not matter anymore, so 
only their coupling with the remaining qubits must be 
refocused. Fig. ED gives an example of such a simplified 
refocusing scheme for five coupled spins. 



1 



J 



1 

2 ± I I 

3 + + 

4 ± ± 

5 ± ± 

FIG. 16 Simplified refocusing scheme for five spins, designed 
such that the coupling of qubits 1-2 with qubits 3 — 5 is 
switched off, i.e. J13, J14, J15, J23, J24 and J25 are inactive 
whereas J12, J34, J35 and J45 are active. 



V 



u 



u 



FIG. 13 Simplification rules for quantum circuits, drawn us- 
ing standard quantum gate symbols, where time goes from 
left to right, each wire represents a qubit, boxes represent 
simple gates, and solid black dots indicate control terminals. 

Sometimes, a quantum gate may be replaced by an- 
other quantum gate, which is easier to implement. For 
instance, refocusing sequences (section IIII.A.40 can be 
kept simple by examining which couplings really need 
to be refocused. Early on in a pulse sequence, several 



More generally, the relative phases between the entries 
in the unitary matrix describing a quantum gate are irrel- 
evant when the gate acts on a diagonal density matrix. In 
this case, we can for instance implement a CNOT simply 
as X2 Uj(l/2J) Y2 rather than the sequence of Eg. 1311 

At the lowest level, that of pulses and delay times, 
further simplification is possible by taking out adjacent 
pulses which cancel out, such as X and X (an instance of 
the first simplification rule of Fig. H5) l. and by converting 
"difficult" operations to "easy" operations. 

Cancellation of adjacent pulses can be maximized by 
properly choosing the pulse sequences for subsequent 
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quantum gates. For this purpose, it is convenient to have 
a library of equivalent implementations for the most com- 
monly used quantum gates. For example, two equivalent 
decompositions of a CNOT12 gate (with J\ 2 > 0) are 



ZyZ-,x, r, ( -L ) v. 



as in Eq. |^ and 



Z 1 Z 2 X 2 Uj[ — ) Y 2 . 



(37) 



(38) 



Similarly, two equivalent implementations of the 
hadamard gate on qubit 2 are 



(39) 



It is thus advantageous to convert as many X and Y 
rotations as possible into Z rotations, using identities 
similar to Eq. 1281 for example 



XY = XYXX = ZX . 



(44) 



and 



A key point in pulse sequence simplification of any 
kind is that the simplification process must itself be ef- 
ficient. For example, suppose an algorithm acts on five 
qubits with initial state 1 00000) and outputs the final 
state (|01000) + |01100)/V2. The overall result of the 
algorithm is thus that qubit 2 is flipped and that qubit 
3 is placed in an equal superposition of |0) and |1). This 
net transformation can obviously be obtained immedi- 
ately by the sequence X 2 Y^. However, the effort needed 
to compute the overall input-output transformation gen- 
erally increases exponentially with the problem size, so 
such extreme simplifications are not practical. 



Y 2 Xi 



(40) 



Thus, if we need to perform a hadamard operation on 
qubit 2 followed by a CNOT12 gate, it is best to choose the 
decompositions of Ecis.l3"7landl4Tjl such that the resulting 
pulse sequence, 



Z t Z, X, F , [ i- I Y, Y, X. 



2 J 



simplifies to 



Z\ Z 2 X 2 Uj 



Xi 



(41) 



(42) 



An example of a set of operations which is easy to per- 
form is the rotations about z. While the implementation 
of z rotations in the form of three RF pulses (Eq. I28|) 
takes more work than a rotation about x or y, rotations 
about z need in fact not be executed at all, provided the 
coupling Hamiltonian is of the form 1\P Z , as in Eq.[2J In 
this case, z rotations commute with free evolution under 
the system Hamiltonian, so we can interchange the order 
of z rotations and time intervals of free evolution. Using 
equalities such as 



ZY = XYXY = XZ, 



(43) 



we can also move z rotations across x and y rotations, 
and gather all z rotations at the end or the beginning of 
a pulse sequence. At the end, z rotations do not affect 
the outcome of measurements in the usual |0), |1) "com- 
putational" basis. Similarly, z rotations at the start of 
a pulse sequence have no effect on the usually diagonal 
initial state. In either case, Z rotations do then not re- 
quire any physical pulses and are in a sense "for free" and 
perfectly executed. Indeed, Z rotations simply define the 
reference frame for x and y and can be implemented by 
changing the phase of the reference frame throughout the 
pulse sequence. 



6. Time-optimal pulse sequences 

Next to the widely used but rather naive set of pulse 
sequence simplification rules of the previous subsection, 
there exist powerful mathematical techniques for de- 
terming the minimum time needed to implement a quan- 
tum gate, using a given system and control Hamiltonian, 
as well as for finding time-optimal pulse sequences l)56l) . 
These methods build on earlier optimization procedures 
for mapping an initial operator onto a final operator via 
unitary transformations l)57t l58f) . as in coherence or po- 
larization transfer experiments, common tasks in NMR 
spectroscopy. 

The pulse sequence optimization technique expresses 
pulse sequence design as a geometric problem in the space 
of all possible unitary transformations. The goal is to 
find the shortest path between the identity transforma- 
tion, /, and the point in the space corresponding to the 
desired quantum gate, U, while travelling only in direc- 
tions allowed by the given system and control Hamilto- 
nian. Let us call K the set of all unitaries k that can 
be produced using the control Hamiltonian only. Next 
we assume that the terms in the control Hamiltonian 
are much stronger than the system Hamiltonian (as we 
shall see in section 1111. B. 21 this assumption is valid in 
NMR only when using so-called hard, high-power pulses). 
Then, starting from /, any point in K can be reached in 
a negligibly short time, and similarly, U can be reached 
in no time from any point in the coset KU, defined by 
{kU\k G K}. Evolution under the system Hamiltonian 
for a finite amount of time is required to reach the coset 
KU starting from K. Finding a time-optimal sequence 
for U thus comes down to finding the shortest path from 
K to KU, allowed by the system Hamiltonian. 

Such optimization problems have been extensively 
studied in mathematics l)59l) . and have been solved ex- 
plicitly for elementary quantum gates on two coupled 
spins l)56tl and a three-spin chain with nearest-neighbour 
couplings l)6Cltl . For example, a sequence was found 
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for producing the trilinear propagator exp(—i2Trlll^I^) 
from the system Hamiltonian h2TrJ(l}l^+I^I^) in a time 
\/3/2J, the shortest possible time l|60l) . This propagator 
is the starting point for useful quantum gates such as the 
doubly-controlled not or toffoli gate. The standard 
quantum circuit approach, in comparison, would yield 
a sequence of duration 3/2 J (it uses only one coupling 
at a time while refocusing the other coupling), and the 
common NMR pulse sequence has duration 1/J. 

Clearly, the time needed to find a time-optimal pulse 
sequence increases exponentially with the number of 
qubits, n, involved in the transformation, since the uni- 
tary matrices involved are of size 2" x 2™. Therefore, the 
main use of the techniques presented here lies in find- 
ing efficient pulse sequences for building blocks acting on 
only a few qubits at a time, which can then be incorpo- 
rated in more complex sequences acting on many qubits 
by adding appropriate refocusing pulses to remove the 
couplings with the remaining qubits. While the exam- 
ples given here are for the typical NMR system and con- 
trol Hamiltonian, the approach is completely general and 
may be useful for other qubit systems too. 



tated by a pulse of length t pw is roughly speaking of or- 
der 1/tpw. Yet, since the qubit response to an RF field is 
not linear (it is sinusoidal in u>it pw ), the exact frequency 
response cannot be computed using Fourier theory. 

For a constant amplitude (rectangular) pulse, the uni- 
tary transformation as a function of the detuning Auj is 
easy to derive analytically from Eas. 1161 and 1171 Alter- 
natively, we can exponentiate the Hamiltonian of Eq. 1151 
to get U directly. An example of a qubit response to a 
rectangular pulse is shown in Fig. 1171 
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B. Experimental Limitations 

Many years of experience have taught NMR spectro- 
scopists that while the ideal control techniques described 
above are theoretically attractive, they neglect important 
experimental artifacts and undesired Hamiltonian terms 
which must be addressed in any actual implementation. 
First, a pulse intended to selectively rotate one spin will 
to some extent also affect the other spins. Second, the 
coupling terms 2TrJijP z I^ cannot be switched off in NMR. 
During time intervals of free evolution under the system 
Hamiltonian, the effect of these coupling terms can easily 
be removed using refocusing techniques fsection lill.ATj) . 
so long as the single-qubit rotations are perfect and in- 
stantaneous. However, during RF pulses of finite du- 
ration, the coupling terms also distort the single-qubit 
rotations. In addition to these two limitations arising 
from the NMR system and control Hamiltonian, a num- 
ber of instrumental imperfections cause additional devi- 
ations from the intended transformations. 



1. Cross-talk 

Throughout the discussion of single- and two-qubit 
gates, we have assumed that we can selectively address 
each qubit. Experimentally, qubit selectivity could be 
accomplished if the qubits are well-separated in space or, 
as in NMR, in frequency. In practice, there will usually 
be some cross-talk, which causes an RF pulse applied on 
resonance with one qubit to slightly rotate another qubit, 
or shift its phase. Cross-talk effects are even more com- 
plex when two or more pulse are applied simultaneously. 

The frequency bandwidth over which qubits are ro- 



FIG. 17 Simulation of the spin response to a 1 ms constant 
amplitude RF pulse as a function of the frequency offset Auj 
between luq and ui r f. The spin starts off in |0) (along +z in 
the Bloch sphere) and wi/27r = 500 Hz is chosen such that 
the rotation angle amounts to 180° for an on-resonance pulse. 

It is evident from Fig. El thai short rectangular pulses 
(known as "hard" pulses) excite spins over a very wide 
frequency range. The frequency selectivity of a pulse can 
of course be increased by increasing t pw while lowering 
B\ accordingly (thus creating what is known as a "soft" 
pulse), but decoherence effects become more severe as 
the pulses get longer. Fortunately, as we will see in sec- 
tions I1V.AI and II V.HI the use of shaped and composite 
pulses can dramatically improve the frequency selectivity 
of the RF excitation. 

Even if a pulse is designed not to produce any net x 
or y rotations of spins outside a specified frequency win- 
dow, the presence of RF irradiation during the pulse still 
causes a shift A.lu* bs in the precession frequency of spins i 
at frequencies well outside the excitation frequency win- 
dow l)6lf) . As a result, each spin accumulates a spurious 
phase shift during RF pulses applied to spins at nearby 
frequencies. 

This effect is related to the Bloch-Siegert shift men- 
tioned in section III. B. II and is known as the transient 
generalized Bloch-Siegert shift in the NMR community. 
It is related to the AC Stark effect in atomic physics. At 
a deeper level, the acquired phase can be understood as 
an instance of Berry's phase l|62T) : the spin describes a 
closed trajectory on the surface of the Bloch sphere and 
thus returns to its initial position, but it acquires a phase 
shift proportional to the area enclosed by its trajectory. 
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The frequency shift is given by 



Alu 



BS 



2(uj - uj rf ) 



(45) 



(provided lu\ <ti \wo — kV/|), where o>o/27r is the original 
Larmor frequency (in the absence of the RF field). In 
typical NMR experiments, the frequency shifts can easily 
reach several hundred Hz in magnitude. We see from 
Eq. 1451 that the Larmor frequency shifts up if lu > uj r f 
and shifts down if loq < Lu r f. 

Fortunately, the resulting phase shifts can be easily 
computed in advance for each possible spin-pulse combi- 
nation, if all the frequency separations, pulse amplitude 
profiles and pulse lengths are known. The unintended 
phase shifts R z {0) can then be compensated for during 
the execution of a pulse sequence by inserting appropri- 
ate R z {— 6), which can be executed at no cost, as we saw 
in section lTlLA.5l 

Cross-talk effects are aggravated during simultaneous 
pulses, applied to two or more spins with nearby fre- 
quencies ujq and Wq (say Wq < u>q). The pulse at Wg then 
temporarily shifts the frequency of spin 2 to + Alubs- 
As a result, the pulse on spin 2, if applied at loq, will be 
off- resonance by an amount —Alubs- Analogously, the 
pulse at Wq is now off the resonance of spin 1 by Alubs- 
The resulting rotations of the spins deviate significantly 
from the intended rotations. 

The detrimental effect of the Bloch Siegert shifts dur- 
ing simultaneous pulses is illustrated in Fig. which 
shows the simulated inversion profile for a spin subject 
to two simultaneous 180° pulses separated by 3273 Hz. 
The centers of the inverted regions have shifted away 
from the intended frequencies and the inversion is in- 
complete, which can be seen most clearly from the sub- 
stantial residual x — ^-magnetization (> 30%) over the 
whole region intended to be inverted. Note also that since 
the frequencies of the applied pulses are off the spin reso- 
nance frequencies, complete inversion cannot be achieved 
no matter what tip angle is chosen (see section llLB.2|) . 

In practice, simultaneous soft pulses at nearby frequen- 
cies have been avoided in NMR (63) or the poor quality 
of the spin rotations was accepted. Pushed by the strin- 
gent requirements of quantum computation, several tech- 
niques have meanwhile been invented to generate accu- 
rate simultaneous rotations of spins at nearby frequencies 
(sections IIV.A.2I and IIV.B.2|) . 



2. Coupled evolution 

The spin-spin couplings in a molecule are essential for 
the implementation of two-qubit gates (section IlII. A. 3|l . 
but they cannot be turned off and are thus also active 
during the RF pulses, which are intended to be just 
single-qubit transformations. Unless loi is much stronger 
than the coupling strength, the interactions strongly af- 
fect the intended nutation. For couplings of the form 
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FIG. 18 Simulation of the spin response to two simultaneous 
pulses with carrier frequencies at Hz and 3273 Hz (vertical 
dashed lines) away from the spin resonance frequency, with 
a calibrated pulse length of 2650^ts (as for an ideal 180°). 
The amplitude profile of the pulses is Hermite shaped (Sec- 
tion IIV.A> in order to obtain a smooth spin response. For 
ideal inversion, the solid line should be —1 at the two fre- 
quencies, and the dashed line should be zero. 



JI\V Z , the effect is similar to the off-resonance effects il- 
lustrated in Fig. 13 the coupling to another spin shifts 
the spin frequency to uo/2ir ± J/2, so a pulse sent at 
u>o/2ir hits the spin off-resonance by T-^/2. 

In practice, J-coupling terms can only be neglected for 
short, high-power pulses used in heteronuclear spin sys- 
tems: typically J < 300 Hz while u>\ is up to » 50 kHz. 
For low-power pulses, often used in homonuclear spin sys- 
tems, uj\ can be of the same order as J and coupling ef- 
fects become prominent. The coupling terms also lead 
to additional complications when two qubits are pulsed 
simultaneously. In general, the qubits become partially 
entangled UK). 

As was the case for cross-talk, NMR spectroscopists 
have developed special shaped and composite pulses to 
compensate for coupling effects during RF pulses while 
performing spin-selective rotations. In recent years, the 
use of such pulses has been extended and perfected 
for quantum computing experiments (sections IIV.AI 
andUVjg. 



3. Instrumental errors 

A number of experimental imperfections lead to er- 
rors in the quantum gates. In NMR, the most common 
imperfections are inhomogeneities in the static and RF 
magnetic field, pulse length calibration errors, frequency 
offsets, and pulse timing and phase imperfections. 

The static field Bq in modern NMR magnets can be 
made homogeneous over the sample volume (a cylinder 5 
mm in diameter and 1.5 cm long) to better than 1 part 
in 10 9 . This amazing homogeneity is obtained by metic- 
ulously adjusting the current through a set of so-called 
"shim" coils, which compensate for the inhomogeneities 
produced by the large solenoid. At ojq — 500 • 2-7T MHz, 
linewidths of 0.5 Hz can thus be obtained, corresponding 
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to a dephasing time constant T 2 * (see Section IV. A. 2\ of 
l/(27r0.5) = 0.32 s. In the course of long pulse sequences 
(of order 0.1 — 1 s), even the tiny remaining inhomogene- 
ity would therefore have a large effect, so its effect must 
be reversed using refocusing sequences f section llII.A.4|) 

The RF field homogeneity is typically very poor, due 
to constraints on the geometry of the RF coils: the en- 
velope of Rabi oscillations (section IV. A. 1(1 often decays 
by as much as 5% per 90° rotation, corresponding to a 
quality factor of only ss 5. In sequences containing only 
a few pulses, this is not problematic, but in multiple- 
pulse experiments, the RF field inhomogeneity is often 
the dominant source of errors and signal loss. 

Imperfect pulse length calibration has an effect similar 
to B\ inhomogeneity: the qubit rotation angle is differ- 
ent than was intended. Only the correlation time for the 
error is different. Miscalibrations are constant through- 
out an experiment, whereas the RF field experienced by 
any given molecule changes on the timescale of diffusion 
through the sample volume. 

Frequency offsets occur in different contexts. In tradi- 
tional NMR experiments, the Larmor frequencies are of- 
ten not known in advance. RF pulses are then expected 
to rotate the spins over a wide range of frequencies, quite 
the opposite case of quantum computing, where the Lar- 
mor frequencies are precisely known and rotations should 
be spin-selective. However, we have seen earlier that P Z P Z 
coupling terms act as a frequency offset of one spin, which 
depends on the state of the other spin. Qubit-selective 
rotations of qubit i thus require a uniform rotation over 
a range ± J2 j7 ti \Jij\/ 2 - 

Various approaches have been developed to reduce the 
sensitivity of RF pulses and pulse sequences to these 
instrumental errors, sometimes in combination with so- 
lutions to cross-talk and coupling artifacts. These ad- 
vanced techniques are the subject of the next section. 



IV. ADVANCED PULSE TECHNIQUES 

The accuracy of quantum gates that can be achieved 
using the simple pulse techniques of the previous sec- 
tion is unsatisfactory when applied to multi-spin systems, 
where the given NMR system and control Hamiltonian 
lead to undesired cross-talk and coupling effects. In addi- 
tion, the available instrumentation can only imperfectly 
approximate ideal pulse amplitudes, timings, and phases, 
for realistic sample geometries and coil configurations, 
and any real molecule includes additional Hamiltonian 
terms such as couplings to the environment, which are 
undesired. Nevertheless, extremely precise control can be 
achieved despite these imperfections, and this is accom- 
plished using the art of shaped pulses, composite pulses 
and average Hamiltonian theory, the subject of this sec- 
ond major section of this review. 

These advanced techniques are based on the assump- 
tion that errors are, at least on some accessible timescale, 
systematic, rather than random. This assumption clearly 



holds for the terms in the ideal NMR Hamiltonian of 
Eqs. [21 an d [23 an d applies also to most instrumental 
errors. Then, by using the special properties of evolution 
in unitary groups, such as the SU(2 n ) which describes 
the space of operators acting on n qubits, the systematic 
errors can in principle be canceled out. 



A. Shaped pulses 

The amplitude and phase profile of RF pulses can be 
specially tailored in order to ease the cross-talk and cou- 
pling effects discussed in sections IlII . B . 1 1 and IIII . B~2l In 
practice, the pulse is divided in a few tens to many hun- 
dreds of discrete time slices; to achieve an arbitrarily 
shaped pulse, it suffices to control the amplitude and 
phase of the slices separately. Furthermore, multiple 
shaped pulses applied at various frequencies can be com- 
bined into a single pulse shape, since a linear vector sum 
of pulse slices also results in a valid pulse. Here, we con- 
sider simple amplitude and phase shaped pulses. 



1. Amplitude profiles 

The frequency selectivity of RF pulses can be much 
improved compared to standard rectangular pulses with 
sharp edges, by using pulse shapes which smoothly mod- 
ulate the pulse amplitude with time. Such pulses are 
typically specially designed to excite or invert spins over 
a limited frequency region, while minimizing x and y ro- 
tations for spins outside this region 

Furthermore, specialized pulse shapes exist which min- 
imize the effect of couplings during the pulses. Such self- 
refocusing pulses lf65h take a spin over a complicated tra- 
jectory in the Bloch sphere, in such a way that the net 
effect of couplings between the selected and non-selected 
spins is reduced (Fig. IT9")l . It is as if those couplings 
are only in part or even not at all active during the pulse 
(couplings between pairs of non-selected spins will still be 
fully active but their effect can be removed using stan- 
dard refocusing techniques IIII. A. 4|l . As a general rule, it 
is relatively easy to make 180° pulses self-refocusing, but 
much harder to do so for 90° pulses. 

The self-refocusing behavior of certain shaped pulses 
can be intuitively understood to some degree. Never- 
theless, many actual pulse shapes have been the result 
of numerical optimizations. Often, the pulse shape is 
expressed in a basis of several functions, for instance a 
Fourier series (165 



(t) = \A + J2 



, 2tt 
A n cos ( n 1 

trim 



B n sin ( n 1 

tnin 



(46) 

and the weights of the basis functions, A n and B n , are 
optimized using numerical routines such as simulated an- 
nealing. 
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FIG. 19 Trajectory on the Bloch sphere of a qubit initially in 
|0), when a so-called iburpI pulse \§5) is applied, of duration 
1 ms and u>\ = 3342 Hz, with a frequency offset (analogous to 
I\I Z coupling) of 0, 100 and 200 Hz. This pulse is intended 
to rotate the qubit from |0) {+z) to |1) (— z). We see that 
the effect of the frequency offset is largely removed by the 
specially designed pulse shape; all three trajectories terminate 
near —z. 



Comparison of the performance of various pulse shapes 
is facilitated by computing the corresponding spin re- 
sponses. This is most easily done by concatenating the 
unitary operators of each time slice of the shaped pulse, 
as the Hamiltonian is time-independent within each time 
slice. Fig. |20| presents the amplitude profile and pulse re- 
sponse for three standard pulse shapes of equal duration, 
illustrating that different pulse shapes produce strikingly 
different spin response profiles. 

Properties relevant for choosing a pulse shape include: 

• frequency selectivity: product of excitation band- 
width and pulse length (lower is more selective), 

• transition range: the width of the transition region 
between the selected and non-selected frequency re- 
gion, 

• power: the peak power required for a given pulse 
length and tip angle (lower is less demanding), 

• self-refocusing behavior: degree to which the J cou- 
pling between the selected spin and other spins are 
refocused (the signature for self-refocusing behav- 
ior is a flat top in the excitation profile) , 

• robustness to experimental imperfections such as 
pulse length errors, 

• universality: whether the pulse performs the in- 
tended rotation for arbitrary input states or only 
for specific input states. 

Table [H] summarizes these properties for a selection 
of widely used pulse shapes. Only universal pulses (also 
known as general-rotation pulses) are included in the ta- 
ble, since quantum computations must work for any in- 
put state. 

Obviously, no single pulse shape optimizes for all prop- 
erties simultaneously, so pulse shape design consists of 
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FIG. 20 (Left) Time profile and (Right) frequency response, 
displaying the z and x — y components of the Bloch vector 
after a pulse when the Bloch vector is along +£ before the 
pulse for three relevant pulse shapes. 
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UBURP 90 
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REBURP180 
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TABLE II Properties of relevant pulse shapes. The Gaus- 
sian |6r?) and Hermite l|67fl shapes are described by analytical 
functions and have been identified early on. The burp family 
of pulses ijBBt) resulted from numerical optimization routines. 
Continued work in this area has produced a large number of 
additional pulse shapes, such as the AV 90 i68fl . 



finding the optimal trade-off for the desired applica- 
tion. For quantum computing experiments, we can select 
molecules with large chemical shifts, so sharp transition 
regions are not so important. Furthermore, the probe 
and spectrometer can deal with relatively high powers. 
The crucial parameters are the self-refocusing behavior, 
the selectivity (short, selective pulses minimize decoher- 
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ence) and to some extent the robustness. 

It is also possible to start from a desired frequency re- 
sponse, and invert the transformation to find the pulse 
shape that produces this response. Again, given the non- 
linear nature of the response, the inverse transformation 
is not given by a Fourier transform, but it can neverthe- 
less be computed directly |o9l) . 

Even self-refocusing shaped pulses do generally not re- 
move the coupling terms completely. Furthermore, when 
two spins are pulsed simultaneously with self-refocusing 
pulses, the refocusing effects are often destroyed |(34|). In 
both cases, the remaining coupled evolution that takes 
place during the pulses must be reversed at an earlier 
and/or later stage in the pulse sequence. 

If we could decompose the evolution during an actual 
pulse into an idealized, instantaneous X or Y rotation 
with no coupling present, followed and/or preceded by a 
time interval of free evolution, we could compensate for 
the coupling effects simply by adjusting the appropriate 
time intervals of free evolution in between the pulses (sec- 
tion [QJ^D- However, Tt r f and Ttj do not commute, so 
such a decomposition is not possible. 

Nevertheless, the coupled evolution can still be un- 
wound to first order iflol fl6j) . when a time interval of 
reverse evolution both before and after the pulse is used: 



,+iHj t/H -i{H rf +H.j) t pul /h +iHj r/h 



,/n 



where r is chosen such that the approximations are as 
good as possible according to some distance or fidelity 
measure (see section lV.CI) . The optimal r is usually close 
but not equal to t pw /2. In comparison, a negative time 
interval only before or after the pulse, 



,+m.j T/h e -i(H rf +Hj) t pw /h , 



is much less effective. 



, /h +iHj t/H 

' 1 



(48) 



2. Phase profiles 

An alternative to amplitude shaping that is often use- 
ful is frequency or phase shaping. One specific phase 
shaping method utilizes fixed, small increments A6 to the 
phase of successive slices of a pulse to achieve an excita- 
tion profile which is centered at a frequency which differs 
from the RF carrier frequency uj r f by A6/At, where At is 
the duration of each time slice. This technique for shift- 
ing the RF frequency is known as phase-ramping l)70l) . 
We can express the effect of phase ramping mathemati- 
cally by replacing Ea. Illlbv 



TL r f/hwi = cos 



sin 



cos 



AS 

Urft + [(pa + -r-t 



( , 

Urft + I <P0 + -j-rt 



Jrf 



A6 
At 



t - 



V (49) 



The use of phase shifts thus permits us to obtain an RF 
field at a different frequency than is generated by the 
signal generator. Furthermore, the displaced frequency 
can be chosen different for every pulse, and can even be 
varied in the course of a pulse. 

A useful application of phase ramping lies in compensa- 
tion for Bloch-Siegert effects during simultaneous pulses, 
where the RF applied at Wq shifts the resonance frequency 
of spin j to u! J + AUbs (Section llII.B.l|l . The rotations 
of both spins can be significantly improved simply by 
shifting the RF excitation frequencies via phase ramp- 
ing such that th ey t rack the shifts of the corresponding 
spin frequencies l)71j) . In this way, the pulses are always 
applied on-resonance with the respective spins. The cal- 
culation of the frequency shift throughout a shaped pulse 
is straightforward and needs to be done only once, at the 
start of a series of experiments. 

Fig. shows the simulated inversion profiles for the 
same conditions as in Fig. 1181 but this time using the 
frequency shift corrected scheme. The inversion profiles 
are clearly much improved and there is very little left- 
over x — y magnetization. Simulations of the inversion 
profiles for a variety of pulse shapes, pulse widths and 
frequency separations, confirm that the same technique 
can be used to correct the frequency offsets caused by 
three or more simultaneous soft pulses at nearby frequen- 
cies. The improvement is particularly pronounced when 
the frequency window of the shaped pulse is two to eight 
times the frequency separation between the pulses UMl- 
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FIG. 21 Similar to Fig. llSl but with frequency shift correction. 



B. Composite pulses 

Another practical method for compensating systematic 
control errors in NMR experiments is the application of a 
sequence of pulses instead of a single pulse. This method 
of composite pulses arises from the observation that con- 
catenation of several pulses can produce more accurate 
rotations than is possible using just a single pulse, due to 
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strategic cancellation of systematic errors and other un- 
wanted systematic effects. Composite pulses work partic- 
ularly well for compensating errors arising from the RF 
field inhomogeneity, frequency offsets, imperfect pulse 
length calibration, and other instrumental artifacts in- 
troduced in Section IIII.B.3I They leverage the ability 
to control one parameter precisely to compensate for the 
inability to control another parameter well. We describe 
two approaches to construction of composite pulses: an 
analytical method, and one employing numerical opti- 
mization. 
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1. Analytical approach 

The three parameters which characterize a hard pulse 
are its frequency offset Aw, phase 0, and area hrfBit pw , 
given by the product of the pulse amplitude B\ and pulse 
duration t pw (Section III.B.1|) . In terms of qubit opera- 
tions, errors in these parameters translate directly into 
errors in the axis n and angle 9 of rotation, such that the 
actual operation applied is not the ideal Rh(0) of Eg . l l2*3"j l . 
but rather, 



R h {6) = cxp 



.f(6,n)-a 



(50) 



where f(9,n) is a function which characterizes the sys- 
tematic error. For example, under and over-rotation er- 
rors caused by pulse amplitude miscalibration or RF field 
inhomogeneity may be described by f(9,n) = 9(1 + e)n, 
while RF phase errors may be described by f(9,n) — 
^[n^cose + n y sine, h v cose — n x sine, h z ], where e is a 
fixed, but unknown parameter. The essence of the com- 
posite pulses technique is that a number of erroneous 
operations are concatenated, varying n and 9, to obtain 
a final operation which is as independent of e as possible. 
This is done without knowing e. 

This technique can be illustrated by considering the 
specific case of linear amplitude errors, in which 



Rn{9) = exp 



i9(l + e)h ■ 8 



(51) 



Let the goal be to obtain R x (ir/2). Using as a measure 
of error the average gate fidelity, defined in Eq. I|113|) , we 
find that F(R x (tt/2),R x (tt/2)) = (2 + cos(e7r/2))/3 « 
1 - vr 2 e 2 /24, so the error is quadratic in e for small e. 
Consider, in contrast, the sequence 



BBl e = R^)R^(2h)R^(it)R x (9) , 



(52) 



about the axis 
= cos" 1 (-9 /4vr) 



where Rip(-) denotes a rotation 
[cos </>, sin </>, 0] , and the choice 4> 
is made. This sequence gives average gate fidelity 
F(R x (tt/2),BB1 90 ) ps 1 - 21vr 6 e 6 /16384, which is much 
better than for the single pulse, even for relatively large 
values of e, as shown in Figure 1221 The operation of the 
BBIqq sequence is illustrated graphically in Fig. 1231 



FIG. 22 Plot of the average gate fidelity between the ideal 
R x (ir/2) and actual unitary transforms R x (n/2) (black line) 
and between the ideal R x (ir/2) and the composite sequence 
BBlgo (blue line), as a function of the fraction of over- 
rotation error e. Note how much higher fidelity the BB1 se- 
quence has (the best possible fidelity is 1) , over a wide range 
of errors. 
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FIG. 23 Illustration of the trajectories of a spin as it trans- 
forms under the BB1 pulse sequence of Eq. l52lL starting 
initially in the |0) state. Three trajectories are shown, in 
which the error is 50% under-rotation (top), zero (middle), 
and 50% over-rotation (bottom). Plotted symbols denote the 
endpoints of each pulse in the sequences. 
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A few comments about this result are in order. This 
result is the best which has been presented in the liter- 
ature to-date |l2 EH 0); currently, no pulse sequence 
which cancels out errors to higher order (for all possible 
initial states) has yet been published. It is also fairly gen- 
eral; BBlg approximates R x (0). Also, while composite 
pulses have been widely studied and employed in the art 
of NMR, this sequence is special in that it is universal 
(also termed fully-compensating or class A): the amount 
of error cancellation is independent of the starting state 
of the spin ifTol l7rf) . Other examples of such universal 
composite pulses are the sequence 

i?6o(180)i?3oo(180)i? 60 (180) , (53) 

which performs a X 2 rotation with compensation for 
pulse length errors, and 

^(385)^ (-320)^(25), (54) 

which performs a Y rotation compensating for off- 
resonance errors and to some extent for pulse length er- 
rors as well. 

Earlier, in the original work which introduced the 
concept of composite pulses into NMR 1)771 : l78i) . only 
limited pulse sequences were known, which only work 
for particular initial states; for example, there is the 
common R x (tt /2)R^ v {'k)R x {'k /2), used to approximate 
R x (w). Figure l2*4l illustrates how this simple sequence re- 
moves the effect of errors in either the rotation angle or 
the rotation axis. 




FIG. 24 Trajectory in the Bloch sphere described by a qubit 
initially in |0), when a composite 180° rotation is applied, con- 
sisting of three imperfect rotations, R x (ji /2)R- y (ji)R x [TT /2). 
(Left) The tip angles are set 0, 5, . . . 20% too short. (Right) 
The pulse is applied off-resonance, with (too — uo r f)/u\ — 
0, 0.05, . . . 0.20. In both cases, the effect of the errors in the 
individual pulses is largely removed by the composite pulse. 

Systematic errors in the coupling strengths can also 
be tackled using composite rotations, in order to obtain 
accurate two-qubit gates. This was shown explicitly for 
the case of Ising couplings l|79|) . 

Similar compensation of slowly-fluctuating errors can 
be achieved during a train of pulses, separated by time 
intervals of free evolution. The simplest instance of such 
a pulse train uses only 180° pulses. Off-resonance ef- 
fects in such pulses can be largely reversed by properly 



choosing the phases of the pulses. For instance, and at 
first sight surprisingly, the errors from off-resonant pulses 
X 2 X 2 roughly add up, while they largely compensate 
each other in X 2 X 2 . This cancellation can be appreci- 
ated via a simple Bloch sphere picture (Fig. I25|) . The 
remaining errors are further reduced for a properly cho- 
sen train of four pulses, X 2 X 2 X 2 X 2 , which performs 
markedly better than X 2 X 2 X 2 X 2 (f80r> . Further reduc- 
tion of the effect of off-resonance errors can be obtained 
by using even longer trains of 180° pulses l|80|) . 




FIG. 25 Trajectory in the Bloch sphere of a qubit initially 
in |0), subject to two consecutive 180° pulses, applied off- 
resonance with (luo — iu r f /lu\ = 0.5. (Left) If the two pulses 
are applied with the same phase (X 2 X 2 ), the qubit is taken 
simply along a circular trajectory through |0), and reaches a 
point near |0); to be precise, the 50% resonance offset makes 
the rotation angle (2 2 + l 1 )/2 2 = ^JTJl larger than 360°. 
(Right) In contrast, if the two pulses are applied with opposite 
phases (X 2 X 2 ), the qubit is left far from |0). 

Evidently, quantum computing sequences are not as 
transparent as just a train of 180° pulses. Surprisingly, 
even throughout a quantum computing sequence, the ef- 
fect of RF inhomogencitics can be removed to a large 
extent (HJ), as illustrated in Fig. [213 After completion 
of a routine involving the equivalent of about 1350 90° 
pulses, the measured amplitude was about 15% of the 
full amplitude. Without removal of the effect of RF in- 
homogeneity, the signal would have been buried in the 
noise very rapidly. 

This level of error cancellation was achieved partly 
due to a judicious choice of the phases of the refocus- 
ing pulses. Nevertheless, a more detailed description and 
understanding of the error operators is needed in order 
to fully exploit the potential for error cancellation in ar- 
bitrary pulse sequences. 

2. Numerical optimization 

The composite pulses we discussed in the previous sub- 
section are designed to compensate for certain types of 
errors (mostly over- or underrotations and frequency off- 
sets), and work even when the exact Larmor frequencies, 
spin-spin coupling strengths and the magnitude of the 
errors are unknown. This is the usual case in NMR spec- 
troscopy. However, in quantum computing experiments, 
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10 15 20 25 30 
# iterations 

FIG. 26 Experimental (error bars) and ideal (circles) ampli- 
tude of | III), as a function of the number of iterations of a 
quantum search algorithm l2t ISlI) , for three qubits, executed 
on 13 CHFBr2. Each iteration contains the equivalent of al- 
most fifty 90° pulses. The dotted lines serve to guide the eye. 
Dashed line: the signal decay for 13 C due to decoherence, 
which represents a lower bound on the decay rate. Solid line: 
the signal strength retained after applying a continuous RF 
pulse of the same cumulative duration per search iteration as 
the pulses in the actual experiment (averaged over 3 spins, 
measured up to 4 iterations and then ex trap olated). Similar 
observations have been reported in Ref. 



detailed knowledge of the system Hamiltonian is usually 
available and can be used to tailor the composite pulses 
to the system specifics, taking the degree of quantum 
control one step further. 

Following the notation of Ref. |83t) . we consider the 
concatenation of a number of rectangular pulses, each 
described by four parameters: the pulse duration r m , a 
constant amplitude lo" 1 , the transmitter frequency uj"^ 
and the initial phase <p m , where m indexes the pulse. 
These parameters may be strongly modulated from one 
pulse to the next 3 Via a numerical optimization proce- 
dure, the values of t™, cj™, lq™^ and <f> m are chosen such 
that the resulting net unitary evolution U ne t is as close 
as possible to the ideal unitary transformation U ideal , ac- 
cording to some fidelity measure (section IV. C|) . 

In practice, the number of time slices in the composite 
pulse is increased starting from one, until a satisfactory 
solution is found. While the fidelity function may have 
many local maxima and finding the global maximum may 
therefore take a long time, suitable algorithms such as 
the Nelder-Mead Simplex algorithm l|84|) often succeed 
in finding a reasonably good solution. Furthermore, the 
optimization routine can incorporate penalties on high 
powers, large frequencies and negative or very long time 
periods, in order to prevent the algorithm from returning 
infeasible solutions. 



Computation of U ne t uses the fact that the Hamilto- 
nian during a fixed-amplitude RF pulse can be made time 
independent by transforming into a reference frame ro- 
tating at the transmitter frequency, as we have seen in 
section lTl.B.2l We will call H" l ot the effective Hamiltonian 
in the frame rotating at lo™^ during segment m. Given 
that uj™f may be different for every segment of the pulse, 
it is most convenient to transform back to a common ref- 
erence frame at the end of every time slice. This can 
be the frame of the raw RF frequency, or the laboratory 
frame of the n-spin system. In the lab frame, the time 
evolution during segment m is described by 



-mi 



(55) 



Since all U m are expressed in the same reference frame, 
we can simply multiply them together to get U ne t — 
Ilm U m , and compare the result directly with U ideal , ex- 
pressed in the laboratory frame as well. 

Two representative examples of composite pulses de- 
signed for spin-selective rotations in homonuclear spin 
systems are given in Fig. |57| The gate fidelity (sec- 
tion ^^^21 obtained with these two pulses is displayed 
in Fig. 123 Naturally, the fidelity is close to unity only 
near the resonance frequencies for which the gate was 
designed to work. 

Composite pulses can thus effectively generate accu- 
rate single- and multiple-qubit Hamiltonians, using de- 
tailed knowledge of the system Hamiltonian, and only 
limited knowledge about the errors. Often, however, 
full knowledge of the system parameters is not available, 
and thus methods beyond composite pulses must be em- 
ployed. 



C. Average-Hamiltonian theory 

The average-Hamiltonian formalism offers a versatile 
framework for understanding how to effectively create or 
remove arbitrary terms in the Hamiltonian by periodic 
perturbations, without requiring full knowledge of the 
system dynamics. The refocusing sequences presented 
in section IIII.A.4I and more general multiple-pulse se- 
quences designed to neutralize the effect of dipole-dipole 
couplings can be explained within this framework. Re- 
duction of full dipole-dipole coupling given by Eq. [S] to 
the simplified forms of Eqs. and [7] can also be under- 
stood with average-Hamiltonian theory. 

Following Ref. 

03), 

we first introduce the Magnus 
expansion and then see how we can modify a time- 
independent Hamiltonian via a time-dependent pertur- 
bation. We use two concrete examples to illustrate the 
concepts. 



Jumps in the transmitter frequency can be conveniently realized 
with phase-ramping techniques; as discussed in Section llV.A.21 
this is done by phase shifting the raw RF excitation in fixed 
increments per time so a different RF frequency is obtained. 



1. The Magnus expansion 

The essence of average-Hamiltonian techniques is that 
the evolution U (t) under a time-dependent Hamiltonian 
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FIG. 27 The ideal RF waveform for two example strongly 
modulated pulses. The solid (dashed) line is the amplitude 
(phase) of the waveform. Details of the pulse parameters, as 
in Eq. 1551 are listed below each waveform. The 6 fis time 
interval with zero RF power before and after the composite 
pulses is needed due to experimental implementation issues. 
The composite pulse in (a) performs a 90° rotation on one 
of the 13 C nuclei of 13 C-labeled Alanine and the pulse in (b) 
performs a simultaneous 180° rotation on two 13 C nuclei of 
13 C-labeled Crotonic acid. Courtesy of D.G. Cory. Repro- 
duced from Ref. 



H(t) can be described by an effective evolution under 
a time- independent average Hamiltonian TL, under two 
conditions (§301 : l43f) : (1) H(t) is periodic and (2) the obser- 
vation is stroboscopic and synchronized with the period 
t c of H(t). 

We can then calculate Tt exactly from 



U(t c ) = exp(-iHt c ) 



(56) 



by diagonalizing U(t c ) and taking the logarithm of the 
resulting eigenvalues Q). 

In practice, it is often more convenient to compute H. 
approximately. Let us assume that Tl(t) is piecewise con- 
stant (the analysis can be easily generalized to the case 
of continuously changing Hamiltonians lj3fjl) ): H(t) = Hk 



for J2o 1 n <t <J2o n, and t c 
U(t c ) = exp(-iH„T n ) . 



En 
Tk, SO 

.exp(-i7Y T ) 



(57) 



a) Three-spin jc/2] pulse 



0.5 





Identity / 


■ / \ 


X 1 

\ 


/ / : \ \ 
// \\ 
'/ , • \,. 


/ 7U/2 / 

/ — ' 







8000 6000 4000 2000 -2000 -4000 -6000 
Spin Frequency (Hz) 
b) Four-spin ti] 12 pulse 



0.5 



/ 












\ / 






\ Identity / 






y 


\ 


\ ~N / 






A 




k 


\ / / 




•A 


J \ 


/ \ 


/ v " 


V / x 










. \ / \ 





6000 4000 2000 -2000 -4000 -6000 -8000-10000 
Spin Frequency (Hz) 

FIG. 28 Gate fidelity of the two example pulses of Fig. 1271 
as the resonance frequency of a test spin is varied. The solid 
(dashed) line is calculated with identity (desired transforma- 
tion) as the intended transformation. The vertical dotted 
lines denote the actual chemical shifts for each spin. Cour- 
tesy of D.G. Cory. Reproduced from Ref. 



Repeated application of the Baker-Campbell-Hausdorff 
relation 
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where 
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(58) 
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■} , (61) 



and so forth. This expansion, called the Magnus expan- 
sion |85|) . forms the basis of average-Hamiltonian theory. 



2. Multiple-pulse decoupling 

Let us consider a pulse sequence of n infinitcsimally 
short pulses Uk separated by time intervals of free 
evolution under the system Hamiltonian 7in, and such 
that U n . . . U2U1 = I (for pulses of finite length, the du- 
ration of the pulses must also be included in the average) . 
The pulses correspond to basis transformations, and we 
can thus describe the system evolution via a sequence of 
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time intervals Tk of free evolution under the Hamiltonian 
Hq (fc), with 

H (0) = Wo, (62) 
W (i) = U^HoUu (63) 
W (2) = U^U^ l H U 2 U x , (64) 

and so forth. Note that the order in which the Uk are 
applied to H.q is reversed and that the Uk themselves are 
reversed as well. If we let t c — X)o Tfc ' then the overall 
transformation U(t c ) is given by 

U(t c ) = exp(-m 0{n) T n ) . . . exp(~m 0(0) T ) . (65) 

We can now use the Magnus expansion of Eq. [2D an d 
Eas. l6Ul6Tl where we replace Hk by Ho(k)> to obtain the 
average Hamiltonian Ho which describes the net time 
evolution during t c . The zeroth order average Hamilto- 
nian is given by 

1 ™ 

n Q) =^12 ■ ■ ■ Uk l H U k ...U x . (66) 

c k=0 

The crux of average Hamiltonian theory is that by prop- 
erly choosing the pulse Uk , we can ensure that H { 0) con- 
tain only the desired terms. 

Sophisticated pulse sequences jg^) can also remove 
undesired contributions from the higher-order terms in 
the expansion, although this is generally harder since 
'HfP , 'H^ ■ ■ ■ contain cross-terms between the various 
Hoik)- The commutators involved in the higher-order 
terms do become smaller for shorter cycle times, though, 
so fast cycles result in better averaging. 

We also point out that pulse sequences which satisfy 

^o(fc) = 7~to(n-k) (67) 

or equivalently 

Uk+i = Ul_ k (68) 

contain no contributions of odd orders to TIq, 

H ( o k) =0 for k = 1,3,5,..., (69) 

and thus perform significantly better than other se- 
quences. 

Let us now illustrate the operation of multiple-pulse 
decoupling via two examples. First, the original multiple- 
pulse sequence for removal of dipole-dipole interactions 
is the wahuha-4 sequence il87l) . 

tX tY 2t Y t X t , (70) 

where the pulses are applied to all qubits involved, r 
stands for free evolution under the system Hamiltonian 
for a duration r, and the unitaries are ordered from right 
to left, as usual. The pulses rotate the Zeeman terms 



in the Hamiltonian from — z to —y, —x, —y, and back 
to — z (see Eqs. I62I64H for a duration r, r, 2r, r and 
r, respectively. The zeroth order average Zeeman term 
is thus oriented along — (x + y + z), with strength scaled 
down by a factor 1 / V3. The dipolar Hamiltonian of Eq.HJ 
goes through the forms [3P Z P Z - P ■ P], [31*13 _ p . p] 
and \3P X P X — P ■ P] for equal durations, and is thus zero 
on average. 

By selectively not pulsing specific qubits, it is also pos- 
sible to reintroduce some of the couplings as desired. In 
Fig- El we already saw explicitly how to do this for P Z P Z 
couplings. 

A second example is an extension of the conventional 
spin-echo sequence fSection IV.A.4|I to three component 
spin-echoes (88). In conventional echo sequences, 180° 
pulses about x or y remove the effect of a Hamiltonian 
of the form c z a z . Now we ask ourselves what sequence 
of pulses would freeze the evolution under a Hamiltonian 
of the form 

H = c x a x + Cyijy + c z a z , (71) 

where c x ,c y ,c z are arbitrary coefficients. We can use 
Eq. 1661 to verify that the sequence 

X 2 t X 2 Y 2 tY 2 Z 2 t Z 2 t , (72) 

or equivalently, after simplification, 

X 2 t Z 2 t X 2 T Z 2 T , (73) 

gives a zeroth order average Hamiltonian H" = 0, and 
thus in effect corresponds to a three-component echo se- 
quence. Another way to show this is to note that 

X 2 H X 2 = +c x a x - CyOy - c z a z , (74) 
Y 2 H Y 2 = -c x o x + CyO y - c z a z , (75) 
Z 2 HZ 2 = -c x a x - CyOy + c z a z . (76) 

Clearly, U + X 2 U X 2 + Y 2 U Y 2 + Z 2 H Z 2 = 0, 
and so the sequence of Eq. 1721 gives, to zeroth order, no 
net evolution. Again, if r is sufficiently short, the higher 
order contributions will be negligible. 

3. Reversing errors due to decoherence 

Can we apply multiple-pulse sequences to reverse the 
effect of interactions of a qubit with degrees of freedom in 
the environment? It is not clear a priori that this is possi- 
ble: standard average-Hamiltonian theory assumes that 
we can manipulate both interacting particles involved, 
for instance via RF pulses. However, we have no control 
of degrees of freedom in the environment. 

Remarkably, it is actually possible to remove the ef- 
fect of unwanted interactions with degrees of freedom in 
the environment, even when applying operations to the 
system only (|89t l90t l91i l92|) , provided the control opera- 
tions are applied faster than the fluctuations. Knowledge 
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about the nature of the interactions can be applied to 
simplify the sequence of decoupling operations, and such 
knowledge can even be experimentally extracted ifflsT ). in 
part using a procedure known as process tomography, de- 
scribed in Section lY.B.2l 

If the fluctuations are faster than the accessible control 
operations, errors can be corrected using quantum error 
correction l|2t l94i I95T) . or they can prevented by encoding 
the qubits in a subspace which is not affected by decoher- 
ence joH I^ . This is discussed further in Section IV.DI 



V. EVALUATION OF QUANTUM CONTROL 

The pulse control methods presented in the last two 
sections can have impressive performance, but this is 
very much contingent on having an accurate model of 
the system under control. A variety of techniques have 
been used in NMR to characterize the system dynam- 
ics and to evaluate the performance of control sequences. 
In this last section, we review some of these techniques, 
beginning with a set of standard experiments to deter- 
mine how quantum a qubit system is, then proceeding to 
tomographic methods for fully characterizing system dy- 
namics, and concluding with fidelity metrics for control, 
and implications these have for scalability to control over 
large systems. 



and the Rabi frequency uj\/2-k is proportional to the am- 
plitude of the control field. Observation of Rabi oscilla- 
tions is usually accepted as a signature of quantum co- 
herent behavior. 

In reality, the envelope of the Rabi oscillation signal 
is always damped, due to decoherence as well as instru- 
mental imperfections; measurement of this decay time is 
useful, and known as a nutation experiment. In NMR, 
the Rabi decay time is often much shorter than the in- 
trinsic phase randomization time constant T 2 , due to 
the inhomogeneity of the RF field driving the Rabi os- 
cillation, across the macroscopic sample. In other sys- 
tems, the Rabi decay time may be longer than T 2 , be- 
cause (1) a long pulse can be seen as a concatenation 
of many 180° pulses, which can have a refocusing effect 
( section |III.A.4(1 . and (2) the qubit is near ±i, where 
phase randomization has no effect, for roughly half the 
time during Rabi oscillations. 

Coherent oscillations driven by a resonant field have 
been observed in NMR and in many atomic systems for 
a long time. Recently, however, observations of such co- 
herent dynamics have been made in other qubit systems, 
including systems made from J osephson junctions l|99|) . 
in molecular vibrational states lllOOt fTTTT^t . and excitons 
in semiconductor quantum dots l|lf)2h . 



2. Coherent oscillations initiated by a kick 



A. Standard experiments 

In NMR spectroscopy as in atomic physics, a num- 
ber of standard experiments serve to test the quantum- 
mechanical behavior of a given system, and to determine 
the extent of its isolation from the environment (see Sec- 
tion III.C|) , in terms of its phase coherence time T2 and 
its energy relaxation time T\, as well as the decay time 
in the rotating frame T\ p . 



1. Coherent oscillations driven by a resonant field 

The dynamics of a single spin, driven resonantly 
by a coherent field, were presented in sections III. B. II 
and lII.B.2l From Eg. 1771 we have that in the ideal case 
the RF field induces transitions from |0) to |1), where a 
qubit initially in |0) will be found in |1) after an RF pulse 
of duration t pw with probability 

Pr[|l>] = sin 2 ( 7 S 1 t p „/2) = sin 2 (u;iW 2 ) , (77) 

The probability initially increases over time, until it 
reaches a maximum Pr[|l)] = I and then decreases again, 
by stimulated emission, a cycle which keeps repeating it- 
self. 

Such oscillations of a two-level quantum system driven 
by a resonant field are known as Rabi oscillations lf98l) . 



A quantum system starting off in a state which is 
not an eigenstate of the (static) system Hamiltonian, 
will precess about the quantization axis of the system 
Hamiltonian, a motion known as Larmor precession (e.g. 
Section IlI.A.lj) . Such a situation could be realized by 
abruptly changing the system Hamiltonian, e.g. by sud- 
dingly applying a strong static field along £ instead of 
along z. Alternatively, and more realistically in NMR, 
Larmor precession can be initiated by suddenly kicking 
the qubit out of the Hamiltonian eigenbasis. For a nu- 
clear spin with Hamiltonian —huj I z (as in Eq. [JJ, this 
is done by applying a 90° RF pulse, causing a transition 
for instance from |0) to (|0) + 11))/ \/2, which initiates the 
time evolution 



\m) = 



e i"ot/2| ) +e^°*/ 2 |l) 

71 



(78) 



as illustrated in Fig. [21 Like Rabi oscillations, the obser- 
vation of Larmor precession is also a signature of quan- 
tum coherent behavior. 

The Larmor precession is also damped, but in contrast 
to the Rabi decay time, the Larmor decay time, termed 
T 2 *, is never longer than T 2 . Usually, T 2 * < T 2 ; in partic- 
ular, for NMR, 
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1 



T* T 2 T> 



(79) 



where T 2 is the dephasing time constant due to static 
magnetic field inhomogeneities or other instrumental im- 
perfections. 
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Larmor oscillations initiated by a kick have been ob- 
served recently in a variety of systems, including those 
driven resonantly mentioned earlier, and in add ition a 
system of charges in coupled quantum dots l)l03l) . The 
oscillations can be observed directly if the measurement 
basis lies in the x — y plane, as is the case in NMR. If the 
measurement takes place along ±z, we must first change 
basis via a second 90° pulse. 



3. Ramsey Interferometry 



The double-pulse experiment 



X t X 



(80) 



where time goes from right to left (as always, for uni- 
tary transformations given in this article) , and r denotes 
a free evolution period, under the evolution operation 
e -iH S y S T/ h_ Th is is known as a Ramsey interference ex- 
periment (jl04f) . Originally, this "method of separated os- 
cillatory fields" was applied to electronic states of molec- 
ular beams traversing through two microwave excitation 
zones. In NMR, two pulses are involved, separated by 
a delay time r. Ramsey interference is most naturally 
described in the rotating frame of the RF. If the qubit 
starts off along z, the first X pulse rotates it to — y. Then 
the qubit precesses about z for a time r. Finally the sec- 
ond X pulse rotates the ±y component of the qubit state 
to =f£. Components along ±x at the end of the interval 
t remain along ±x after the second X pulse. 

If only a single qubit is considered and the RF field 
is exactly on-resonance with the qubit precession, the 
qubit stays in place in the rotating frame during the time 
interval in between the two pulses, and the final state 
does not vary with r. However, if the RF and the qubit 
are detuned in frequency by Aui, both the x and the z 
components of the final state display a beating pattern 
as a function of Aw r, the so-called Ramsey fringes. The 
decay time of the envelope of the Ramsey fringes is T 2 , 
the same as that for Larmor precession. 

For coupled qubits, the beating pattern contains infor- 
mation on the coupling strengths. This fact forms t he ba - 
sis for two-dimensional correlation spectroscopv (j30t Il05l) . 
a widely used range of two-pulse techniques for molecular 
structure determination. 



4. Measurement of T2 

The intrinsic T2 time can be extracted in an experi- 
ment which is based on Larmor or Ramsey experiments. 
Certain imperfections which cause the Ramsey or Lar- 
mor decay time T 2 * to be smaller than T2 can be removed 
by applying refocusing sequences. 

The simplest instance of such a refocusing sequence 
consists of a single 180° pulse applied halfway during the 
time interval of free evolution initiated by an initial 90° 



pulse. The entire sequence is thus 



- X z - X . 



(81) 



A second X pulse should be added at the end if the mea- 
surement takes place in the ±z basis. In multi-spin sys- 
tems, the pulses must be applied selectively to one spin, 
in order to measure the T2 of that spin. 

The X 2 refocusing pulse removes not only simple scalar 
spin-spin couplings, as described in section IIII.A.4I but 
also undoes the effect of spatial variations of the static 
magnetic field along z. Such field inhomogeneities make 
spins in different regions of the sample progressively get 
out of phase with each other during the first time in- 
terval r/2. As a result, their magnetic moments cancel 
each other out and the NMR signal vanishes. Provided 
the magnetic field variations are constant throughout the 
experiment, all the spins get exactly in-phase again (now 
along +y) by the end of the second time interval r/2, be- 
cause of the 180° refocusing pulse. As a result, the signal 
recovers, producing the well-known spin-echo. A gener- 
alization of this technique known as three- component re- 
focusing fSection |lV.C.2|l . undoes effects from any static 
spin Hamiltonian terms. 

The echo signal decays as a function of r, and the decay 
time constant is a measure of Ti. However, terms in the 
Hamiltonian fluctuating on a timescale shorter than r are 
not removed by a single refocusing pulse. Their effect can 
still be removed if the fluctuations are slow compared to 
r/n and a train of n refocusing pulses is applied, each 
preceded and followed by a time interval r/2n of free 
evolution. This so-called Carr-Purccll sequence (106), 
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(82) 



produces a first echo along +y after r/n, a second echo 
along — y after 2r/n, a third along +y after 3r/n, and 
so forth. The magnitude of the echo signal decays ex- 
ponentially throughout this sequence and the echo signal 
left at the end of this sequence decreases exponentially 
as a function of the total time r. To the extent that slow 
fluctuations in the Hamiltonian have been refocused, the 
decay time constant is the intrinsic T 2 . 

As we have seen in Section IIIII of this review, small 
but fixed errors in the pulse amplitude or duration may 
accumulate throughout a multiple-pulse sequence such as 
the Carr-Purcell sequence. However, if the phase of the 
refocusing pulses is shifted by 90° with respect to the 
initial 90° pulse, pulse length errors are compensated on 
even-numbered echoes and are thus not cumulative. In 
this sequence, 



2n 



Y 1 



,Y 2 



— Y — X 

n 2n 



(83) 



known as th e Carr-Purcell-Meiboom-Gill or CPMG se- 
quence l|l07|) . the echoes appear all along —y. Again, the 
decay time constant of the echo signal gives T 2 . 

Since T 2 indicates for how long a qubit can remain 
phase coherent, it is usually called the coherence time, 
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although the terms phase randomization time, dephasing 
time and transverse relaxation time, are also used. In 
NMR, T 2 is also known as the spin-spin relaxation time. 
In any case, T 2 is an important number for evaluating the 
potential of quantum computers, as the ratio of T 2 over 
the typical duration of a quantum logic gate expresses the 
number of operations that can be completed coherently. 

5. Measurement of Ti 

Energy exchange with the environment makes a qubit 
which is out of equilibrium gradually return to thermal 
equilibrium. In thermal equilibrium, the qubit is in a 
statistical mixture of |0) and |1), with probabilities set 
by the temperature and the energy difference between |0) 
and |1). The time constant of this equilibration process, 
Ti, is often called the energy relaxation time, the longi- 
tudinal relaxation time or simply the relaxation time. In 
NMR, T\ is often termed the spin-lattice relaxation time. 

Two standard experiments for measuring T± are inver- 
sion recovery and saturation recovery. The sequence for 
the inversion recovery experiment is 

X t X 2 . (84) 

The 180° pulse inverts the |0) and |1) probabilities, then 
during time r, relaxation takes place, and finally a 90° 
read-out pulse is applied if necessary (i.e. when the mea- 
surement basis is in the x — y plane). In saturation recov- 
ery, a strong RF field is applied for a long enough time 
such that it saturates the qubit transition and equalizes 
the |0) and |1) probabilities. As in inversion recovery, the 
original |0) and |1) populations are altered, and we can 
monitor the populations return to their equilibrium value 
as a function of r. The time constant of this equilibration 
process is T\. 

Note that both the inversion recovery and saturation 
recovery experiments bring the qubit out of equilibrium, 
but to a state which has no coherence. As a result, phase 
randomization does not affect the equilibration process 
— we measure purely the effect of energy exchange with 
the bath. In contrast, Ramsey and spin-echo experi- 
ments for measuring T 2 * an d T 2 pick up contributions 
from phase randomization both without and with energy 
exchange with the bath. If energy exchange dominates 
phase randomization, the measured T 2 is 2Ti. 

The relevance of T\ is twofold. First, it sets an upper 
bound for T 2 and second, it tell us how much time we have 
to perform a measurement in the {|0), |1)} basis. Phase 
randomization does not change the |0) and |1) probabil- 
ities, so T 2 is irrelevant during such a measurement. In 
many cases, Ti 3> T 2 , in which case we have more time 
to measure than to perform coherent operations. 

6. Measurement of Tip 

A third decay time useful to characterize the degree of 
isolation between a qubit and the environment is T\ p . 



This time constant can be measured via a technique 
called spin-locking, where the spin is first rotated into 
the x — y plane, say by a Y pulse, and next continuous 
irradiation is applied, phase shifted by 90° with respect 
to the pulse, so it is aligned with the spin state (along 
the x axis): 

R x (continuous) Y . (85) 

The continuous irradiation along x locks the spin to the 
x axis, in the following sense. Whenever the spin starts 
to diverge from the x axis due to interactions with the 
environment, the RF field rotates it to the opposite side 
of the x axis within a time tt/uji. Provided that the spin 
is still moving in the same direction after this time, it 
will thus return to the x axis (note that spin-locking thus 
also inhibits evolution due to J-couplings and moderate 
frequency offsets). Only if the spin moves in the opposite 
direction after tt/oji, it continues to depart from the x 
axis due to the rotation by the RF field. So the amplitude 
along x decays, and the decay time constant is termed 
Tip, known in NMR as the spin-lattice relaxation time in 
the rotating frame. 

We see thus that, whereas T 2 is governed by low- 
frequency fluctuations in the environment and T\ de- 
pends on fluctuations at the Larmor frequency, the decay 
during spin-locking arises from fluctuations at the Rabi 
frequency used during spin-locking. The spin-locking ex- 
periment thus gives additional information on the spec- 
tral density of the interactions with the environment. 



B. Measurement of quantum states and gates 

The standard experiments presented in the previous 
section give only partial information on the system dy- 
namics. Here we show that in fact the full relaxation 
superoperator can be determined systematically by a pro- 
cedure known as process tomography, which builds upon 
state tomography, as described below. 



1. Quantum state tomography 

The density matrix p completely describes our knowl- 
edge of the state of a system. Measurement of the den- 
sity matrix is therefore extremely helpful when testing or 
claiming the preparation of specific quantum states. 

One-time measurement of each of n qubits, in a given 
basis of 2™ states |m), gives only very little information on 
p. All that can be inferred from a measurement outcome 
to is that Pr[|m)] ^ 0. 

Repeated measurements of n qubits, each time pre- 
pared in the same state and measured in the same basis, 
reveals the probability distribution for the measurement 
basis states, 

Pr[|m)] = (m\p\m) = Tr(p|m)(ra|) = Tr(pM) , (86) 
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where M is an observable or measurement operator. If 
we repeatedly measure each qubit in the {|0), |1)} basis, 
we thus obtain all the diagonal en tries of p. pa. 

Quantum state tomography flit llOSt Il09|) is a method 
which allows all the elements of the density matrix p to 
be determined. This method consists of repeating the 
measurement of the same state in various measurement 
bases, until all the elements of p can be determined, by 
solving a set of linear equations. In practice, it is often 
more convenient to first rotate the qubits via a unitary 
transformation and then perform the measurement in a 
fixed basis. This is equivalent to measuring in different 
basis, since 



Tr [p{UMU r )] = Tr [(U f pU)M] 



(87) 



Specifically, we can expand the density matrix of a 
single qubit p as 



poo|0)(0|+poi|0)(l|+pio|l)(0|+pn|l)(l| 



Poo Poi 
pio pn 



Measurements of a single qubit in the {|0), |1)} basis give 
us poo an d pn = 1 — Poo- However, after changing basis 
via a 90° rotation about x, transforming p to XpX\ 
the measurement gives access to Im(pio) = — Im(poi)' 
Similarly, measurement after transformation by Y reveals 
Re(pio) = Re(poi)- Thus, by measuring the qubit state 
first directly, then measuring the same state again after 
an X read-out pulse, and then again after a Y read-out 
pulse, we can reconstruct p completely. 

Similarly, for n qubits, we can expand p as 



E E Pw'NXil > 



(89) 



i=0 j=0 



and choose a set of basis changes which gives access to 
all 4™ — 1 degrees of freedom in p. 

However, it is much easier to find a suitable set of 
basis changes if we use the Pauli expansion of p instead 
of Eq. The Pauli expansion for a single-qubit state is 



P = CoOo + C\<J\ + C 2 (72 + C3CT3 , 



(90) 



where Co = 1 for normalization, and we use <jq = 1/2, 
&i = cr x /2, (T2 = 03 = <J z /2. Measurement in the 

computational basis, described by the observables ao±a 3 , 
gives us Pr(|0» = (c + c 3 )/2, and Pr(|l» = (c - c 3 )/2 
so we can extract C3. Since 



XpX^ = Co (To -I- C\G\ — C3CT2 + C2O3 
YpY f = C CTo + C3CT1 + C 2 (T2 - C1CT3 



(91) 
(92) 



we indeed obtain (co ± Ci)j2 after applying X and (cq =p 
ci)/2 after using Y . 

For n qubits, Eq. |OT| generalizes to 



3 3 

EE- 

i=0 j=0 



3 

E 

k=0 



.k Oi 



(93) 



where C00...0 = 1- Measurement in the computational 
basis is described by observables of the form 

{a ± a 3 ) (g> (an ± a 3 ) ® . . . {a ± <r 3 ) , (94) 
and returns the probabilities 



E 



± 



i,j,...fc£{0,3} 

For example, for two qubits, these are 

PrflOO)) = (c 00 + c 03 + c 30 + c 33 )/4 

Pr(|01)) = (c 00 - c 03 + c 30 - c 33 )/4 

PrfllO)) = (coo + c 03 - c 30 - c 33 )/4 

Pr(|ll)) = (coo - c 03 - c 30 + c 33 )/4. 



(95) 



(96) 
(97) 
(98) 
(99) 



After measurement of the four Pr[|m)], we can solve for 
c 03j c 30i c 33 from this overdetermined set of linear equa- 
tions. Again, we can determine the other Cy...fc by trans- 
formation of the corresponding (Jij...k to an observable, 
for instance 

A 1 y 2 (cr + cr 2 ) <g> (ct + vi)x\Y% = (cr + 0-3) ® {<r - 03) ■ 

(100) 

We end this discussion of state tomography with three 
additional comments: 

First, in order to obtain all the basis state probabilities 
such as in Eqs. 1961991 we must each time read out all 
the qubits. If it is only possible to read out any one 
single qubit in each experiment, we obtain n bit-wise 
probabilities instead of 2™ basis state probabilities, giving 
spin-spin correlations. The measurement operators are 
then of the form 

2 ,i ~ 1 [er ® o-o ® (<>"0 ± 0-3) ® . . . O cr ] (101) 
and we measure probabilities 

2 ( c o...ooo...o ± C0...030...0) (102) 

It is now no longer possible to rotate arbitrary compo- 
nents of p into observable positions using just single-qubit 
rotations. Two-qubit gates are necessary to obtain all 

Cij ...k- 

Second, the measurement basis need obviously not be 
the computational basis. In NMR experiments, for in- 
stance, the single-qubit measurement operator can be 
written as — ia\ — 02- For two coupled spins, the mea- 
surement operators are 

2(-t<7l-£T 2 ) ® (^0±cr 3 ) (103) 

{<J ±a 3 ) ® 2(-icr 1 - a 2 ) , (104) 

and so forth. Since NMR experiments are normally per- 
formed on a large ensemble of molecules, the expectation 
value of the observables can be read out by acquiring 
a single spectrum. The four operators in Eq. 11041 cor- 
respond to the four lines in the spectrum of a two-spin 
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system (two doublets) . Phase sensitive detection permits 
us to separately record the real and imaginary component 
of each spectral line and distinguish <j\ and 02 contribu- 
tions. 

Third, errors in the gates used for changing basis dur- 
ing state tomography, lead to a measured density matrix 
which differs from the actual state of the system. If the 
errors are known and can be modeled accurately, they 
can be incorporated in the state tomography procedure 
and the actual state can nevertheless be determined ac- 
curately. 

Quantum state tomography has been experimentally 
implemented in many atomic syste ms, notably the early 
work mapping out photon st ates ijllOft and vibrational 
cat states of trapped atoms Recently, it has be- 

come a common tool used to evaluat e NM R states ifiTI: 
llOSft . states of optical photon qub its l)ll2j) . and even vi- 
brational states of molecules ill 13ft. 



2. Quantum process tomography 

Now that we know how to experimentally determine 
the state of a quantum system, it is only a short step 
to the characterization of a quantum process, such as 
a quantum logic gate, communication channel, storage 
device and so forth. In general, let us consider a quantum 
mechanical black box whose input may be an arbitrary 
quantum state, and whose output is the result of the 
internal dynamics of the black box, as well as interactions 
to the outside world. Then can we ascertain the transfer 
function of this black box? 

The answer is yes JH1 El E3 [TT3). The outline 
of the procedure is to first determine the output state 
of the black box for a set of input states which form 
a basis for the system Hilbert space, and then to use 
the fact that quantum mechanics is linear to compute 
the entire transfer function from this finite set of input- 
output pairs. 

An arbitrary quantum state transformation is a linear 
map £, 



£{p) 

Tv{£{p)) 



(105) 



where we can express £{p) in the oper ator-sum represen- 
tation or Kraus representation Q IllSft (an alternative to 
the superoperator formalism widely used in NMR H3)), 



(106) 



The Afc are operators acting on the system alone, yet £ 
completely describes the possible state changes of the sys- 
tem, including unitary operations, generalized measure- 
ments and decoherence (for trace-preserving processes, 
J2i A\A k — 1). The expansion of Eg. 11061 is in general 
not unique. In fact, we can always describe £ using a 



fixed set of operators A k which form a basis for the set 
of operators on the state space, so that l|ll4ft 



£(P) = XI X ^A 



vP 



(107) 



where Xvq ^ s a positive Hermitian matrix. Since the A k 
are fixed, £ is completely described by X- I n general, \ 
will contain 16™ — 4™ independent real parameters, where 
n is the number of qubits. 

In order to determine x experimentally, we choose a 
basis of 4™ linearly independent density matrices pj which 
span the system Hilbert space, and determine £(pj) for 
each j. We can then write down a set of linear equations 
of the form of Eg. 117771 where we plug in the measurement 
outcomes £ (pj) and solve for the Xmn- 

The most convenient choice for the pj depends on 
the implementation of the qubits and on the observ- 
ables, as was the case for state tomography. Clearly, 
the effort needed to perform quantum process tomog- 
raphy increases even more rapidly with the number of 
qubits than quantum state tomography. This procedure 
has ex perimenta lly been used only fo r one and two qubit 
NMR i|ll7t ITTgft and optical photon IpSfm^ systems. 

The operation elements A k in the operator sum rep- 
resentation of Eq. 11061 can describe arbitrary quantum 
operations, but among these a select subset are useful 
to identify. For example, when A k is a unitary matrix, 
that corresponds to perfect, closed system Hamiltonian 
evolution. Phase damping (T2) is described by 



A = 



1 

v/7 







% /T^ 



(108) 



and amplitude damping (T±) by 



A 



1 

V7 



Ai 



x/T - 




(109) 



where 7 ~ e~ l l r parameterizes the strength of the damp- 
ing, for time t, in terms of a time constant r. These, and 
other relaxation parameters 10) , can be obtained by pro- 
cess tomography. 

Such results can, in turn, be useful for approximate 
numerical simulation of relaxation and decoherence pro- 
cesses in spin systems. Phase damping and energy relax- 
ation can be simulated in alternation with unitary evo- 
lution under the system and control Hamiltonian, taking 
sufficiently short time slices to obtain good approxima- 
tion of true dynamics. This permits a n spins system to 
be simulated using n4™ steps, compared to 16™ for fully 
general quantum operations. Experimental results have 
shown this method to be predictive of the system dy- 
namics thr oug hout sequences containing hundreds of RF 
pulses ifll l^. 
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C. Fidelity of quantum states and gates 

The methods of the previous section give us full knowl- 
egde of the system state and dynamics, but sometimes a 
more succinct measure for comparing theoretical expec- 
tations with experimental measurements is useful. These 
are given by quantum state and gate fidelities. 



1. Quantum state fidelity 

One elementary goal of quantum control is to create 
some pure state However, suppose the final output 
is instead the pure state |c/>). Does \<j>) represent \if>) with 
high fidelity? 

Classically, the fidelity of two probability distributions 
{p x } and {q x } is given by F(p x ,q x ) = J2 X y/PxQx', when 
they are equal, the fidelity is one. The analogous quan- 
tum measure of fidelity for two pure states \ip) and \<f>) 
is 

F(M,|0)) = KM>>|, (no) 

which is simply the absolute value of the overlap between 
the two states. 

More generally, the output state of a control sequence 
is often described by a density matrix p; this is useful 
because density matrices can describe classical statisti- 
cal mixtures of quantum states, arising from decoherence 
processes, for example. The fidelity between a pure state 
\tp) and a mixed state p is 

nw.rf-vW, (in) 

which reduces to EcL ljllUfl when p — |(/>)(</>|. 

The most general case is the fidelity between two den- 
sity matrices, p and er, which is defined as (0) 

F(a,p)=Tv S J^pV^, (112) 

and despite the apparent asymmetry in this expression, 
it is actually symmetric in p and er, and furthermore, 
reduces properly to Eq. 1 (111( 1 when one density matrix is 
pure. 

Note that in the literature, s omet imes the square of 
Eq. Ijl 10(1 is defined as the fidelity l|l22() ; that departs from 
the usual classical definition for fidelity, but is convenient 
because F(\ip), \4*)) 2 can be interpreted as the probabil- 
ity that a system in the state \<fi) is found to be in the 
state \ip) when measured in the {\ip), \il>±)} basis. Such 
probabilities are meaningful in the accuracy thresholds 
discussed in Section IV. Dl 

Other metrics for comparing two states have been used 
to quantify the relative error between theoretical and ex- 
perimental states, such as the simple two-norm (IsTh and 
other expressions |83j) . These were used because the di- 
agonal elements of the density matrix were suppressed; 
such metrics are inferior to the fidelity measure, which 



should be used when possible, due to its direct connec- 
tion to quantum information measures and fault toler- 
ance theorems. 

It is worthwhile to consider a specific example re- 
lating control precision to state fidelity. Suppose we 
desire \ip) — |1), but obtain \<f>) = Ry(ir + e)|0) = 
-sin | |0)+cosf |1) « -e/2|0) + (l-e 2 /8)|l). The result- 
ing error probability is 1 — |(0|V')| 2 = e 2 /4. This example 
makes the point that that for small rotation angle errors 
e, the gate failure probability goes as e 2 . 



2. Quantum gate fidelity 

A more complex goal of quantum control is to accom- 
plish a desired quantum operation. Perhaps the most 
common scenario is one in which the desired operation U 
is a unitary transform on a single qubit, whereas the ac- 
tual transform accomplished is some quantum operation 
£ (given in the operator sum representation). 

A natural means to evaluate control precision is 
through the average gate fidelity 

F(£,U) = J F{\^)^£{\^){i>\)Ufd^ (113) 
= J (ii\tf£(\ij)(i>\)U\i>)di>, (114) 

where the integral is over the uniform (Haar) measure dtp 
on the Hilbert space of the system. For a single q ubit , 
this formula can be reduced to a simple expression (|X22t 
OH, 

F{£,U)= l ~ + ± Tr(Ua k U^£(a k )), (115) 

fc={l,2,3} 

where a k are the three Pauli matrices. Similar simple 
formu las c an be obtained for higher dimensional sys- 
tems (123). Note that by convention (0), the average 
gate fidelity is defined such that it goes as the square of 
the usual state fidelity; thus, it can be interpreted as a 
probability. 

A more difficult quantity to calculate is the minimum 
gate fidelity, 

F(£ : U)=mm w F(U\i;),£(\i>)(m ; (116) 

the square of this quantity gives the worst-case gate fail- 
ure probability that is relevant for fault-tolerance thresh- 
old theorems. 



D. Evaluating Scalability 

This article has been concerned with the control of 
complex systems composed of multiple distinct physical 
pieces. Given some degree of control over a few such 
pieces, how controllable is a very large quantum system 
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composed from many pieces? Normally, one would ex- 
pect that a system composed of unreliable pieces would 
itself be unreliable, and that the overall probability of 
failure increases rapidly with the number of pieces. Unex- 
pectedly, however, arbitrarily reliable quantum systems 
can be built from unreliable parts as long as certain cri- 
teria are met. 

The main criterion for being able to construct a reliable 
system is that the probability of error per oper ation p be 
below the "accuracy threshold," VL2& llM liH iHl 
Pth- When p < pth is satisfied, a quantum error cor- 
rection circui t (for ins t ance as demonstrated by NMR 
experiments (Il29t Il30t Il3ll) ) can be constructed using 
the unreliable components; this circuit performs compu- 
tations on encoded qubits, such that a net decrease in 
error is achieved even when error correction itself is done 
with the faulty gates. 

The probability of failure per operation p, must of 
course be defined, and this is done in terms of the fidelity 
metrics discussed in the previous section, which incorpo- 
rate decoherence (e.g. T\, T2, gate times, etc.) and con- 
trol imperfections. Thus, for example, p is bounded from 
above by the gate fidelity p < 1 — F{£, U) 2 . 

Remarkably, no reliable resources need be utilized for 
the fault tolerant construction. Through k levels of re- 
cursive application of error correction, the device error p 
can be reduced to p 2 , using physical resources (space, 
time, and energy) which scale as d k for some constant 
d. Thus, a small increase in resources exponentially re- 
duces the overall error. Many assumptions are made in 
obtaining p t h, such as the availability of local, fast, paral- 
lel classical control resources, but the generall y accepte d 
theoretical optimal value oipth is about 1CP 4 (126; 132), 
with optimistic estimates ranging as high as 10~ 3 with 
additional restrictions As we have seen at the end 

of section IY.C.1I this implies that fo r insta nce rotation 
angles must be precise to order ~ VlO~ 4 = 10~ 2 . In 
principle, p t h can be experimentally measured, for exam- 
ple, by implementing a recursive error correction circuit 
and testing its probability of failure, but this has not yet 
been accomplished. 

The fault tolerance threshold pth, and its relative value 
compared with state and gate fidelities give a crisp crite- 
rion for system scalability for specific implementations. 
Modern classical systems are robust mainly because com- 
ponent failures can be controlled; similarly, the future of 
control over quantum systems hinges on our ability to 
evaluate p t h and to build components which fail with 
probability p < p t h- 



VI. DISCUSSION AND CONCLUSIONS 

In this review, we have presented a diverse set of tools 
intended to compensate for undesired or uncontrolled 
terms in the Hamiltonian of coupled qubits, as well as for 
instrumental limitations. These tools are most powerful 
and easiest to design when all the terms in the system 



Hamiltonian commute with each other, and the control 
terms can be much stronger than the system Hamilto- 
nian. The common theme of the control techniques is 
careful tailoring of the amplitude, phase and frequency of 
the time-dependent terms in the Hamiltonian, whether in 
the form of shaped pulses, composite pulses or multiple- 
pulse sequences. We now discuss the effectiveness and 
applicability of these advanced control techniques, which 
points at where they could be used in other quantum 
systems. 

Pulse shaping is particularly attractive because of 
the modular and scalable design approach. Amplitude 
profiles are selected from a library of standard or spe- 
cially designed shapes in order to minimize cross-talk 
(frequency-selective pulses) and coupling effects (self- 
refocusing pulses). Robustness to experimental imper- 
fections can also be considered in the choice of pulse 
shape. Once suitable amplitude profiles have been cho- 
sen, the pulse lengths are set as short as possible while 
maintaining qubit-selectivity. The same amplitude pro- 
files and pulse lengths are then used throughout the pulse 
sequence. 

Remaining cross-talk effects can be further reduced 
at a small cost (quadratic in the number of qubits). 
From the amplitude profiles and pulse lengths, unin- 
tended phase shifts produced by single RF pulses as well 
as off-resonance effects during simultaneous pulses can 
be precomputed, once for every pair of qubits. 

The main disadvantage of the standard pulse shaping 
techniques is that often the coupled evolution during the 
pulses (in particular 90° pulses or simultaneous pulses) 
cannot be completely frozen. The remaining coupled evo- 
lution can be unwound to a large extent during the time- 
intervals before and after the pulse, but such reversal is 
never perfect because the RF terms in the Hamiltonian, 
P x and do not commute with the coupling terms, l\P z - 
Furthermore, shaped pulses must often be quite long in 
order to remain spin-selective, which means that deco- 
herence has more effect. This problem evidently gets 
worse as the Larmor frequencies of the spins approach 
each other. 

Nevertheless, the combination of pulse shaping and 
phase ramping techniques has been very successful in 
practice. It has enabled the implementation of the most 
complex sequences of operations realized to date, acting 
on up to seven nuclear spins. 

Composite pulses have proven to be a versatile tool 
in NMR spectroscopy, mostly for compensating system- 
atic errors such as RF field strength variations and fre- 
quency off-sets. Another useful application of composite 
pulses is the effective creation of unitary operators which 
are otherwise not or not easily accessible. A good exam- 
ple is the composite z rotation, created from a sequence 
of x and y rotations. 

Even so, the use of (hard) composite pulses in NMR 
quantum computing experiments has been limited so far. 
Their main drawback is that in multi-spin homonuclear 
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molecules, single-frequency but high-power and rectan- 
gular pulses will rotate spins in a large frequency win- 
dow, about an axis and over angles which depend on RF 
field strength and the respective resonance off-sets. This 
severely limits straightforward application of hard com- 
posite pulses in homonuclear spin systems. 

Nevertheless, it is in principle possible to take advan- 
tage of the differences in resonance off-sets in order to 
rotate one spin while the other spins undergo no net ro- 
tation. Such effective frequency selectivity despite the 
use of hard pulses was demonstrated in a quantum com- 
putation on a hom onuc lear two-spin system, first using 
single hard pulses l)l34h . and later using composite hard 
pulses Jz2). 

The same idea underlies the operation of composite 
pulses tailored to achieve any rotation of one or more 
spins about independent axis, using detailed knowledge 
of the system Hamiltonian. Furthermore, short, high- 
power pulses can be used, so the effect of decoherence 
is reduced compared to the case of the long, low-power 
shaped pulses. Even more attractive here is the fact 
that all the coupling terms can be effectively frozen, and 
that other types of cross-talk, such as Bloch-Siegert ef- 
fects, are automatically taken care of, sunlike the case of 
shaped pulses. 

The main disadvantage of such strongly modulated 
composite pulses is that the time needed to find near- 
optimal pulse parameters increases exponentially with 
the number of qubits n, as it involves computing unitary 
matrices of size 2™ by 2™. Nevertheless, for small 
numbers of qubits, this technique can be very useful. 

Average-Hamiltonian techniques underlie the oper- 
ation of widely used multiple-pulse refocusing sequences. 
In the context of liquid NMR quantum computing, cou- 
plings are of the form P Z P Z and refocusing sequences con- 
sist simply of a train of 180° pulses. Such refocusing se- 
quences are an essential ingredient of all NMR quantum 
computing experiments involving more than two spins. 

More complex decoupling sequences exist to remove 
the effect of coupling Hamiltonians of a different form, 
as is the case of solid-state NMR and many other qubit 
implementations. Even errors arising from interactions 
with the environment, i.e. decoherence, can be removed 
using multiple-pulse sequences. 

In all cases, the refocusing operations (e.g. the 180° 
pulses) must be fast compared to the fluctuations they 
are intended to cancel, and they must also be repeated 
at a rate faster than the fluctuations. 

Perspective — In NMR quantum computing experi- 
ments on heteronuclear spin systems, where short, high- 
power RF pulses were used, errors in the time evolution 
were usually dominated by various instrumental limita- 
tions. Most experiments on homonuclear systems, in con- 
trast, made use of long, low-power pulses, and were lim- 
ited by cross-talk and coupling effects. As the pulse tech- 
niques for coping with limitations of the Hamiltonian and 



instrumentation became more advanced, the field reached 
the point where errors due to imperfect quantum control 
were smaller than errors caused by decoherence. 

Reaching this point in many-qubit systems must be 
a prime objective for any implementation of quantum 
computers, along with reduction of decoherence itself. 
Quantum information and computation theory offers a 
common language which can facilitate transfer and trans- 
lation of the techniques for coherent control of coupled 
nuclear spins to other fields of physics. Such cross- 
fertilization has alrea dy started, in systems as di verse 
as trapped ions l)l35|) . excit ons in quantum dots <|l36f) 
and Cooper pair boxes l)l37[) . and is likely to accelerate 
the progress towards the elusive goal of complete control 
over quantum systems. 
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